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We investigate the effects of a K-mouflage modification of gravity on the dynamics of clusters of 
galaxies. We extend the description of K-mouflage to situations where the scalar field responsible for 
the modification of gravity is coupled to a perfect fluid with pressure. We describe the coupled system 
at both the background cosmology and cosmological perturbations levels, focusing on cases where 
the pressure emanates from small-scale nonlinear physics. We derive these properties in both the 
Einstein and Jordan frames, as these two frames already differ by a few percents at the background 
level for K-mouflage scenarios, and next compute cluster properties in the Jordan frame that is better 
suited to these observations. Galaxy clusters are not screened by the K-mouflage mechanism and 
therefore feel the modification of gravity in a maximal way. This implies that the halo mass function 
deviates from A-CDM by a factor of order one for masses M > 10 14 K~ 1 Mq. We then consider the 
hydrostatic equilibrium of gases embedded in galaxy clusters and the consequences of K-mouflage 
on the X-ray cluster luminosity, the gas temperature, and the Sunyaev-Zel’dovich effect. We find 
that the cluster temperature function, and more generally number counts, are largely affected by 
K-mouflage, mainly due to the increased cluster abundance in these models. Other scaling relations 
such as the mass-temperature and the temperature-luminosity relations are only modified at the 
percent level due to the constraints on K-mouflage from local Solar System tests. 

PACS numbers: 98.80.-k 


I. INTRODUCTION 

K-mouflage M is one of the four types of screen¬ 
ing mechanisms, together with the chameleon I'i. Q, the 
Damour-Polyakov Q and the Vainshtein Q ones, which 
are compatible with second-order equations of motion 
for single scalar field models [?']. The properties of K- 
mouflage have already been thoroughly studied both at 
the background cosmology |7j and perturbation levels 
[ 1 , ;9|, see also 0 for a more general analysis within an 
“imperfect-fluid” formalism. The small-scale dynamics 
have been studied in 0, and models that satisfy both 
cosmological and Solar System constraints have been de¬ 
vised in 121. 

In this paper, we extend these studies by including flu¬ 
ids with pressure as befitting the description of gases in 
galaxy clusters. We do so for both the background and 
perturbations. We also present the dynamics of the sys¬ 
tem in both the Einstein frame (used in previous works) 
and the Jordan frame, and discuss the relations between 
both frames. Because the properties of gases, such as the 
X-ray luminosity or the Sunyaev-Zel’dovich effect fl3i | . 
or the wavelength of atomic emission or absorption lines, 
are tied to the frame in which atomic physics is described 
without any modification due to the scalar field, we work 
in the Jordan frame to describe clusters of galaxies. We 
use the fact that galaxy clusters are not screened by the 
K-mouflage mechanism and that their number would be 
increased as compared with A-CDM in this scenario. As 
a result, clusters appear as a useful testing ground for 
K-mouflage and its effects on the growth of structure. 

We take into account the modification of the hydro¬ 


static equilibrium in K-mouflage models, together with 
the change of the matter density profiles, which we find 
to become slightly more compact. This allows us to inves¬ 
tigate the X-ray luminosity as a function of the gas tem¬ 
perature. The deviation from A-CDM is at the percent 
level and is set by the tests of gravity in the Solar System, 
which strongly constrain the coupling constant that de¬ 
fines these models. In a similar fashion, the temperature- 
mass relation is affected at the same level. As particu¬ 
lar examples, we focus on a simple “cubic” K-mouflage 
model (that agrees with cosmological constraints) and 
on an “arctan” model which satisfies cosmological con¬ 
straints as well as all Solar System tests^ including the 
advance of the perihelion of the Moon [12| |. The latter 
gives slightly amplified results as compared to the for¬ 
mer, but both remain at the percent level. The only ob¬ 
servable which deviates significantly from A-CDM is the 
cluster temperature function (or more generally, number 
counts) due to the increased abundance of clusters for 
masses M > 10 14 h~ 1 M & . 

The paper is arranged as follows. In section un we de¬ 
fine the K-mouflage models and detail some of their prop¬ 
erties for a fluid with pressure coupled to K-nrouflage at 
the background and perturbation levels, working in the 
Einstein frame. In section Hill we reformulate the dynam¬ 
ics in terms of the Jordan-frame quantities. In section llVl 
we present our numerical results for the background and 
the growth of structure for two K-mouflage models: the 
cubic and arctan models. In sectionEl we turn to galaxy 
clusters and their properties focusing on the physics of 
the gas embedded in the clusters. In section ED we dis¬ 
cuss in details the similarity and differences between the 
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K-mouflage scenarios and other modified-gravity theo¬ 
ries. We conclude in section rvm 

A derivation of the equations of motion in the Ein¬ 
stein frame is given in appendix s while details on 
the Einstein-Jordan connection can be found in the ap- 
pendixQjJ We discuss unitarity constraints in appendix IU1 


II. DEFINITION OF K-MOUFLAGE MODELS 

A. Jordan-frame and Einstein-frame metrics 


We consider scalar-field models where the action has 
the form 00 


S = 




R + C v {ip) 


the Friedmann equations are modified). In this paper 
we compute the properties of astrophysical objects such 
as clusters of galaxies, including their temperature and 
X-ray luminosity, and it is more convenient to work in 
the Jordan frame. Then, radiative processes, such as 
Bremsstrahlung, are given by the standard results and 
do not vary with time or space. Moreover, matter den¬ 
sity is conserved. This simplifies the analysis, as the only 
difference from a A-CDM scenario will be a change of 
gravity laws, which can be explicitly derived from the 
action HU). 

The conformal transformation © actually means that 
the line elements are transformed as ds 2 = A 2 ds 2 . Using 
conformal time r and comoving coordinates x, this local 
change of distance can be absorbed in the scale factor for 
the background universe, as 

ds 2 = a 2 (—dr 2 + dx 2 ), ds 2 = a 2 (—dr 2 + dx 2 ), (3) 


+ /d^£-(«&W), (1) 

which involves two metrics, the Jordan-frame metric Vl 
with determinant g , and the Einstein-frame metric g 
with determinant g. The matter Lagrangian density, £ m , 
where ipm are various matter fields, is given in the Jordan 
frame, where it takes the usual form without explicit cou¬ 
pling to the scalar field (although one could add explicit 
couplings to build more complex models). The gravita¬ 
tional sector is described by the usual Einstein-Hilbert 
action, but in terms of the Einstein-frame metric g^ 

and the associated reduced Planck mass Mpi = 1 /\/ 8 nQ. 
The Lagrangian density C v {p) of the scalar field is also 
given in the Einstein frame. 

Throughout this paper, we denote Einstein-frame 
quantities with a tilde, to distinguish them from their 
Jordan-frame counterparts (when they are not identical). 
We choose this notation, which is thejrpposite to the one 
used in our previous papers 0 0 [lj, [lj| where we mostly 
worked in the Einstein frame, as here we mostly work in 
the Jordan frame. 

If the two metrics were identical, this model would be a 
simple quintessence scenario [lj, ll5| , with an additional 
scalar field to the usual matter and radiation components 
but with standard electrodynamics and gravity (General 
Relativity). In this paper, we consider modified-gravity 
models where the two metrics are related by the confor¬ 
mal transformation [16] 

9nu=A 2 (<p)g tlv . (2) 

This gives rise to an explicit coupling between matter 
and the scalar field. In the Einstein frame we recover 
General Relativity (e.g., the Friedmann equations) but 
the equations of motion of matter are non-standard (e.g., 
the continuity equation shows a source term and matter 
density is not conserved). In the Jordan frame the equa¬ 
tions of motion of matter take the usual form (e.g., mat¬ 
ter density is conserved) but gravity is modified (e.g., 


with 

a = Ad , r = f, x = x. (4) 

[Throughout this paper, we denote with an overbar mean 
background quantities, such as A = A{(p)] However, 
physical time t and distances r, with ds 2 = —dt 2 + dr 2 , 
are changed as 

dt = Adt, r = ax = Ar. (5) 

In particular, the cosmic times t and t are not the same 
in both frames 


B. K-mouflage kinetic function 

In this paper, we consider K-nrouflage models 000, 
which correspond to cases where the scalar-field La¬ 
grangian has a non-standard kinetic term, 

£ (p (ip) = M 4 K(x) with x= - 7 ^ 4 ^ (6) 

Throughout this paper, (V^) is the covariant deriva¬ 
tive associated with the metric gp v ( g (hence x = 
A~ 2 x, but we work with \ in the following). Here, Af 4 
is an energy scale that is of the order of the current dark- 
energy density, (i.e., set by the cosmological constant), to 
recover the late-time accelerated expansion of the Uni¬ 
verse. Thus, the canonical cosmological behavior, with 
a cosmological constant pa = M 4 , is recovered at late 
time in the weak-y] limit if we have 

X->0: K(x)~- 1 +(7) 

where the dots stand for higher-order terms, the zeroth- 
order factor —1 corresponding to the late-time cosmo¬ 
logical constant M 4 . The normalization of the first two 
terms in Eq. J7]) defines the normalizations of the constant 
AA 4 and of the field ip, hence it does not entail any loss of 
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generality (within this class of models). We only consider 
models that satisfy this low-x expansion in this article, 
and where \ —■ y 00 for t —> 0 and x —)• 0 for t —> 00 . 

Well-behaved K-nrouflage scenarios have K' > 0, 
where we denote K' = dK/d'x. and W±(y ) = 

yK'(±y 2 / 2) must be monotonically increasing functions 
up to +00 over y > 0. This ensures that the cosmological 
dynamics are well defined up to arbitrarily high redshift, 
where the matter density becomes increasingly large, and 
that small-scale static solutions exist for any matter den¬ 
sity profile [11 . Moreover, there are no ghosts around 
the cosmological background nor small-scale instabilities 
0 - 


We must point out that the kinetic functions K (x) 
that we use for numerical computations and illustra¬ 
tive purposes in this paper are defined by fully nonlin¬ 
ear expressions, namely Eqs. CD and m below, and as 
such go beyond the low-x expansion 0- As explained 
above, the latter expansion is very general and holds 
for well-behaved models, where K' > 0 for all x and 
W±(y) = yK' (±y 2 /2) are monotonically increasing func¬ 
tions of y. The expansion 0 would only be violated if 
K' diverges at low x, e.g. K{x) = — 1 + x 3 / 4 + •••, but 
we do not consider such singular cases here. 

Then, it happens that at at low redshifts, in the dark- 
energy era, x [with its normalization defined by the first 
two coefficients in the expansion 0] is small on cosmo¬ 
logical scales, which implies K' cs 1. This holds both 
for the homogeneous background and for the cosmologi¬ 
cal large-scale structures. This property is related to the 
fact that at low redshifts, in the dark-energy era, we re¬ 
quire the cosmological evolution to remain close to the 
A-CDM behavior. From the expressions m , we can see 
that this implies xK' -C K (to recover a dark-energy 
equation of state pd e — —Pde) whence x€l. In fact, at 
low z we have the scaling x ~ /3 2 , where j3 is the cou¬ 
pling strength defined in Eq.0 below, so that x ~ 0.01 
as we take /? = 0.1. We shall check this behavior in Fig. 0] 
below. 


We shall also check in Sec. IVBl and Fig. [0] below that 
this also applies to clusters of galaxies at low redshifts, 
which are not screened by the nonlinearities of the scalar- 
field Lagrangian, in spite of their large mass. This would 
not be the case for a coupling /3 0.1, but this would 

violate some Solar System and cosmological constraints 
and we do not consider such models here. 

Nevertheless, the nonlinearities of the kinetic function 
K (x) come into play at high redshift and are taken into 
account in our computations, using the explicit nonlinear 
examples ED and ED- This ensures in particular that 
the dark-energy density becomes subdominant at high 
2 and that we recover the Einstein-de Sitter cosmology 
in the early matter era 0]. Moreover, the background 
solution can be shown to be stable and is a tracker solu¬ 
tion 0. The nonlinearities on the far negative semiaxis, 
—X 1, also play a critical role to ensure that Solar 

System tests of gravity are satisfied by the K-mouflage 
model, but we do not consider this regime in this paper. 


Although K-mouflage theories involve high-order 
derivative interactions, they do not suffer from quantum- 
mechanical problems such unitarity violation in their in¬ 
teraction with matter [TtI [Ts[ , as explained in the ap¬ 
pendix [0 


C. K-mouflage coupling function 


The coupling function A(ip) has the low-tp expansion 

A(<p) = 1 + + ..., ( 8 ) 

M pi 

where the dots stand for higher-order terms. The nor¬ 
malization of the first term does not entail any loss of 
generality and only corresponds to a normalization of 
coordinates. At early times, t —> 0, we have <p —> 0 and 
9 [iv —>• More generally, we define the coupling /? as 


/%) = Mpi 


din A 

d tp 


(9) 


It is constant for exponential coupling functions, A(ip) = 
exp[/3(^/Mpi]. Without loss of generality, we take /3 > 0 
(which simply defines the sign of the scalar field ip). 

Cosmological and Solar System constraints imply /3 < 
0.1, see Ref. 0. Moreover, we have the scaling 
\/3tp/Mpi\ ~ /3 2 <C 1, see Ref.0, as we shall check 
in Fig. [4] below (see also Ref. [9[). Therefore, in re¬ 
alistic models, we have \A — 1| < 0.1 and the higher- 
order terms in the expansion 0 only have a small quan¬ 
titative impact. We shall also check in Fig. [13] be¬ 
low that the fluctuations of the scalar field are small, 
\<p — (p\ <C |<p|, so that the coupling function A(ip) re¬ 
mains dominated by the low-order terms of the expan¬ 
sion @ in clusters of galaxies (and at smaller scales). 
This can be readily understood from the fact that real¬ 
istic models should have a fifth force that is not greater 
than the standard Newtonian force. This typically im¬ 
plies |<L4/A| < |Tn|, where H/n is the Newtonian poten¬ 
tial, whence \(38tp/Mpi\ < 10 -5 . 


D. Equations of motion in the Einstein frame 

Observable effects, such as lensing or two point correla¬ 
tions that can be measured, are independent of the choice 
of frame, so that we can work in either the Einstein or 
the Jordan frame. As explained in the introduction, for 
our purposes the Jordan frame is more convenient and 
more transparent. Indeed, in this frame both the matter 
and radiation components obey their usual equations of 
motion, e.g. the matter energy-momentum tensor sat¬ 
isfies = 0 so that the matter density obeys the 

usual conservation equation. Moreover, particle masses 
and atomic emission or absorption lines do not evolve 
with the cosmic time (whereas they do in the Einstein 
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frame). Then, the only effect of the scalar field is to 
change the gravitational sector, that is, the Friedmann 
equations that determine the background cosmological 
expansion rate and the relation between the metric grav¬ 
itational potentials and the matter density fluctuations 
(i.e., it leads to modified Poisson equations that can be 
interpreted as a fifth force). 

Therefore, in this article we work in the Jordan frame 
and compute observable effects, in particular the prop¬ 
erties of clusters of galaxies, in this frame. However, 
to simplify the derivation of the equations of motion, 
it is convenient to first derive the Friedmann equations 
and the equations that govern the growth of cosmolog¬ 
ical large-scale structures in the Einstein frame, where 
gravity takes the standard form. In a second step, we 
will use these results to obtain the equations of motion 
in the Jordan frame through a change of variables, in 
Sec. EH Afterwards, all our computations will remain in 
the Jordan frame. 

Thus, we describe in the appendix E the derivation 
of the equations of motion of the scalar field and of the 
matter component in the Einstein frame, for a cosmolog¬ 
ical fluid with a nonzero pressure. In this section we only 
give the main results, which will be needed to obtain the 
equations of motion in the Jordan frame in Sec. Hill 

We consider three components of the energy den¬ 
sity of the Universe, a matter fluid with nonzero pres¬ 
sure, radiation, and the scalar field. The Einstein-frame 
and Jordan-frame matter energy-momentum tensors are 
given by 


T = 

/I v — 


-2 SS m 

y/^gSg^ 


T — 

fllS - 


-2 SS m 
V~9 ’ 


( 10 ) 


where we dropped the subscript “m”. The conformal 
transformation © gives 


T lu , = A-*T liV , T^ = A~ 4 f^ T tlv = A~ 6 T^, (11) 


small-scale (sub-horizon) limit, k/aH 1, where spatial 
gradients dominate over time derivatives and quasi-static 
approximations apply. Thus, we focus on the regime de¬ 
fined by 

v 2 <1, T n « 1, | ~ c? ~ § N « 1, > 1. (13) 

p aH 

As the background level, the matter background pres¬ 
sure is zero, p = 0, and the Einstein-frame Friedmann 
equations read as 

3MpjJJ 2 = p + p^ + pip, (14) 

d JJ 

= p + P(i) + P(i) + Pip + Pip, (15) 

where p, /5( r ) and pq) = /5( r ) /3 are the background matter 
and radiation densities and pressure, which evolve as 


dp 

dt 


_dln A 


dp(r) 

dt 


dJJp(r), 


(16) 


while p v and p v are the background scalar-field density 
and pressure, given by 

P v = - K), % = M 4 K, (17) 


with 


~ = _J_ 

X 2A4 4 ^ dt ) 


(18) 


(Throughout this paper, we consider a flat universe with 
zero background curvature.) The background scalar field 
obeys the Klein-Gordon equation 


d_ 

dt 


-3 

a —-K 
dt 


, -din .A 


= —a p- 


d(p 


(19) 


where we use g^ (g^ u ) to raise indices in T (T), and the 
relation g^ = A~ 2 g tiv . In particular, the Einstein-frame 
and Jordan-franre densities and pressures are related by 

P = A 4 p, p = A 4 p. (12) 

We work in the non-relativistic limit, v 2 <C 1 where 
v is the mean fluid peculiar velocity, and in the weak 
field regime, 4 /n 1, where 4% is the Newtonian grav¬ 

itational potential. Moreover, assuming the usual Cold 
Dark Matter (CDM) scenario for the dark matter, the 
matter pressure p is negligible on cosmological scales and 
it arises from the small-scale nonlinear processes, such 
as the collapse of gas clouds which generate shocks, or 
the virialization of of dark matter halos (which generate 
an effective pressure through the velocity dispersion of 
the matter particles). Then, p ~ pc 2 , where c s is the 
speed of sound or the velocity dispersion, and c 2 ~ Tn 
because it is generated by the gravitational collapse (as 
in hydrostatic equilibrium where pressure gradients bal¬ 
ance the gravitational force). In addition, we consider the 


For the matter perturbations, the continuity and Euler 
equations write as 

dp , - din A , 

_^ + V .(pv) + 3Hp = p^ r , (20) 

dv , _ dlnA\ ~ Vp 

— + (v-V)v+ ( % -I--j— j v = -V(4/ N + lnA)-—, 

(2i) 

where the Einstein-frame Newtonian potential is given 
by the Poisson equation 

= <22) 

In the small-scale (quasi-static) limit, the Klein-Gordon 
equation for the scalar field becomes 

1 _ H In 4 

-V • (Vp K') = ——5p. (23) 

a z dp 
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Here we also used the fact that the fluctuations of <p can 
be neglected in the factor K' : so that the Klein-Gordon 
equation can actually be linearized in the scalar field, 
while keeping the matter density fluctuations nonlinear. 
See Ref. Q for a detailed discussion and an explicit com¬ 
putation of the matter power spectrum, up to one-loop 
order, that includes up to the cubic term in <p in Ea. (E3ll . 
which is checked to make no quantitative difference for 
cosmological large-scale structures of cluster sizes and be¬ 
yond. This corresponds to the fact that clusters are not 
screened by the nonlinear K-mouflage mechanism, which 
comes into play at much smaller scales and higher densi¬ 
ties, as in the Solar System. 


III. EQUATIONS OF MOTION IN THE 
JORDAN FRAME 

We now derive the equations of motion of the scalar 
field and of the matter component in the Jordan frame. 
To do so, we use the results obtained in Sec. Ill Dl in the 
Einstein frame, and express these equations in terms of 
Jordan-frame variables. 


A. Background dynamics 


The Jordan-frame metric g^ is related to the Einstein- 
frame metric g „„ by the conformal transformation (f 2 |) . 
As seen in Eqs.j3|)-([5]), this leads to a rescaling of the 
scale factor and of physical time and distance, while the 
conformal coordinates are unchanged. The Hubble ex¬ 
pansion rates, H = dlna/df and H = dlna/dt, are also 
different and related by 


H{ l + e 2 ) H 

A ~ A( l-e 2 )’ 


(24) 


where £2 (t) was defined in Eci. (IA45l) and verifies 
dlnA/dt = egJJ, and we introduced its Jordan-frame 
counterpart, 


£2 {t) 


d In A £ 2-62 

d In a 1 + £2 1 — £2 


(25) 


The Einstein-frame and Jordan-frame densities and 
pressures are related as in Eq. G2D. so that the Friedman 
equation m gives 


in the Jordan frame, as compared with the usual Gen¬ 
eral Relativity result, because the gravitational Einstein- 
Hilbert action is defined in terms of the auxiliary metric 
g^. Substituting for g^ this effectively corresponds to 
a change of the Einstein-Hilbert action. At the back¬ 
ground level, this simply introduces the time-dependent 
functions A(t) and £2 (t) in Eg. (1251) . 

We can also write Eq. (OBI) in the standard form (albeit 
with a time-dependent reduced Planck mass), as 

3Mp { H 2 = p + /5( r ) + pdej (28) 

by defining the dark-energy component as the energy 
density that is “missing” in the Friedmann equation to 
match the Hubble rate, after we sum over the other mat¬ 
ter and radiation components. This yields 

2 g 2 _ £ 2 

Pde = Pip + 7~j v2 (.P Pv’) ' (29) 

(i - £ 2 j 

This interpretation corresponds to the case where mea¬ 
surements of the Hubble rate and of the matter and ra¬ 
diation densities are performed in the Jordan frame and 
the remaining part, which explains the accelerated ex¬ 
pansion, is ascribed to the dark-energy component, [as 
in the usual A-CDM case, where the background dark 
energy is also measured from the missing energy density 
that is required to account for H(z)\. This is a natural 
configuration, as a(t) and H(t) are obtained from redshift 
measurements of standard candles, which assumes that 
atomic absorption and emission lines are the same at dis¬ 
tant redshifts as in the laboratory. By definition, this is 
the case in the Jordan frame but not in the Einstein frame 
(where particle masses are actually time dependent). On 
the other hand, these standard candles must not depend 
on the local gravity, because in the Jordan frame New¬ 
ton’s constant becomes time dependent, so that these 
candles are no longer standard (i.e., similar to those at 
z = 0). This rules out supernovae (which involve both 
local gravity and electrodynamics, within the star) but 
allows one to use geometric candles such as baryon acous¬ 
tic oscillations [Hi [20| or the Alcock-Paczynski test [2l| . 

Using Eq. (usd, the matter and radiation densities now 
evolve as 

t = - 3 "A TT = - 4H ^y < 30 > 

Thus, we recover the usual conservation equations for 
matter and radiation in the Jordan frame, whence 


3Mp X H 2 — (1 — £ 2 ) 2 (p + P ( r ) + Pip) ■, (26) 

where we introduced the Jordan-frame reduced Planck 
mass, 

M 2 l {t) = M 2 l /A(i) 2 . (27) 

Thus, in the Jordan frame, Newton’s constant, Q — 
l/8nMp l: varies with time, as Q(t) = QA 2 oc A 2 . Equa¬ 
tion (EO) shows how the Friedmann equation is modified 


p( t ) = - 1 > P(r)(t) = ^-, with o 0 = l, (31) 

a J a 4 

where po are the mean Jordan-frame energy densities to¬ 
day, at z = 0, and we normalized the Jordan-frame scale 
factor by ao = 1 . 

From Ea. IAllI) . with p = 0, and Ea. (l29l) . the Jordan- 
frame dark-energy density evolves as 

dpde 
df 


3 H (pde T Pde) > 


(32) 
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where we defined the Jordan-frame dark-energy pressure 
as 


where we introduced the Jordan-frame Newtonian poten¬ 
tial, given by 


Pde 


= P<p + 


£2 

1-62 


(P(r) +Pv) + 



2 de 2 \ 
1 — e 2 d In a J 


P + P(r) + Pep 

3(1 -e 2 ) 2 


(33) 


On the other hand, from Eq. © the Klein-Gordon equa¬ 
tion reads as 


d_ 

df 


A-V^K' 

d t 


= — a 3 p 


ddn A 

dtp 


This can be integrated as 


fir = -A\ 
d t 



PoP(t') 

Mpi 


(34) 


(35) 


because the integration constant must vanish to recover 
a realistic early-time cosmology |7). 

Finally, we define the Jordan-frame cosmological pa¬ 
rameters as 


On = —, fi (r) = -^-, Oie = —, (36) 

Pcrit Pcrit Pcrit 

where /9 cr ; t = 3Mp X H 2 = A _4 (l — e 2 ) _2 /5 cr it is the Jordan- 
frame critical density. This gives 

fim = (1 - e 2 ) 2 D m , n (r) = (l-e 2 ) 2 H (r) , (37) 


and 


^de — &ip + (2e 2 — e 2 ) (f2 m + fi(r))- (38) 


-4v 2 ^n = 77 VjT S Pt whence T N = ^ N . (42) 

CL Zildpj 

The last equality follows from Eas. (B^ll and (E71) . There¬ 
fore, in the Jordan frame the two metric potentials are no 
longer equal, but their sum remains equal to 24 , N . This 
is related to the fact that photons do not feel the effect of 
the fifth force, see also Eq. (1A71) . Therefore, weak lensing 
statistics show the same dependence on the matter den¬ 
sity fluctuations as in GR and the impact of the modified 
gravity only arises through the different evolution of the 
density field and the time-dependent Newton constant, 
see also Sec lV HI below. 

The Klein-Gordon equation (l23l) reads as 



pA 

K'Mpi 


Sp, 


and, up to first order over Sip, we obtain 


SA 

A 


P 

Mp\ A 


5tp, 


1 „2 SA 

—V 2 —— = 

a 2 A 


Mpj/f' 


Sp. 


(43) 


(44) 


This also gives 

=Isifr 5 '’ “ d <45) 

where e\ is defined by 


ei (t) 


2/3 2 

K” 


(46) 


We can check that +flde = f2 m + ^(r) + D V = 1. 

The effective dark-energy equation of state in the Jordan 
frame is simply defined as 

Hide — Pde/ Pde- (39) 


B. Perturbations 

The dynamics of large-scale perturbations in the Jor¬ 
dan frame are also obtained from the equations derived 
in the Einstein frame in Sec. GU2 

In the Einstein frame, the Newtonian gauge metric 
(IA16I) reads as ds 2 = a 2 [— (1 + 2^N)dr 2 + (1 — 24 , N)dx 2 ], 
where we used Ea. (IA31l) . In the Jordan frame, we write 

ds 2 = a 2 [—(1 + 2<f>)dr 2 + (1 - 2T)dx 2 ]. (40) 

Then, using ds 2 = A 2 ds 2 and a = Aa, we obtain, up to 
first order in SA = A — A, 

-r T T T JA 

$ = 4'n + ^-, ^ = (41) 


and we recover the same factor as in the first relation 

(1X451) . 

The continuity equation (12011 and the Euler equation 

m become 

|^ + V • (pv) + 3U P = 0 (47) 

and 

<9v V v 

— + (v- V)v + Hv = -V4>--. (48) 

OT p 

Therefore, in contrast with the Einstein frame, in the 
Jordan frame the continuity and Euler equations take the 
same form as in A-CDM, and the coupling to the scalar 
field tp only gives rise to the modified Poisson equation 
(1451) . in terms of the formation of large-scale structures. 
There is no longer a non-conservation term in the con¬ 
tinuity equation nor an additional friction term in the 
Euler equation. However, in contrast with the Einstein 
frame and the A-CDM cosmology, the two gravitational 
potentials $ and T that enter the Newtonian gauge met¬ 
ric are now different. 
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C. Formation of large-scale structures 


D. Spherical collapse dynamics 


Introducing the Jordan-frame matter density contrast, 
6 = Sp/p, (49) 

the continuity equation (l47l) also writes as 

+ V • [(1 + <5)v] = 0. (50) 

This is the same equation as (IA40I) and we have 

J = <5 = <5. (51) 

Indeed the matter density contrasts in the Einstein 
frame, whether we consider the density p or the “con¬ 
served” density p of Ea. (IA13l) (in case of zero pressure), 
and in the Jordan frame, are equal within the approxi¬ 
mations described in the previous sections. 

On perturbative scales, we again set the pressure term 
to zero and we introduce the two-component vector 


(£) “ (-(V-v)/(da/dt)) ' (52) 

Equations (IA42I) and (IA43I) become 

cflna ~ = J dkldka <5D(ki+k 2 -k)d(ki,k 2 ) 

x V’ 2 (ki)'i/'i(k 2 ), (53) 


As can be derived from Ea. (l48l) . on large scales where 
the pressure is negligible, the particle trajectories r(f) 
read as 


d 2 r 

di 2 


1 d 2 a 

ad? r 


-v r $, 


(56) 


where r = ax is the physical coordinate and V r = V/a 
the gradient operator in physical coordinates. To study 
the spherical collapse before shell crossing, it is conve¬ 
nient to label each shell by its Lagrangian radius q or 
enclosed mass Iff, and to introduce its normalized radius 
y{t) by 


m wMl 5 = (Sk) ■ v(t =0)=L 

(57) 

In particular, the matter density contrast within radius 
r(t) reads as 


1 + 8(<r)=y(t) . 


(58) 


Since the Poisson equation (l46l) is only modified by the 
time-dependent prefactor 1 + ei(f) and the time depen¬ 
dence of Newton’s constant, for a spherical system the 
gravitational force is still set by the total mass within 
radius r, 


dj '2 

d In a 


-fl m (l + 


1 d H\ 
JP~At ) 


4>2 = 


dd» . GSM 
— = (1 + ei)——, 
dr r z 


(59) 


fdkidk 2 <5_D(ki+k 2 -k)/3(ki, k 2 )^ 2 (ki)^ 2 (k 2 ). (54) 


We recover the same form as for the A-CDM cosmol¬ 
ogy, except for the factor e± that corresponds to a time- 
dependent amplification of Newton’s gravity, from the 
modified Poisson equation (H5l) . 

On large scales or at early times, we can again linearize 
the equations of motion and the evolution equation (IA46I) 
for the linear modes becomes 


d 2 D 
d(lna) 2 


„ 1 AH 

+ lH~dt 


AD 
din a 


—H m (l + £i )D — 0. 


(55) 

Again, as compared with the Einstein-frame Eq. (IA46I) we 
find that the coefficient e 2 has disappeared and the only 
difference from the A-CDM case is the time-dependent 
amplification of the gravitational term by (1 + ei). 

As in Galileon models, but in contrast with f(R ), Dila- 
ton and Symmetron models, the linear growing modes 
remain scale independent as in the A-CDM cosmology. 
This is because we did not include a potential V(<p) 
in the scalar-field Lagrangian © and the field is mass¬ 
less. Thus, the amplification of gravity extends up to 
the Hubble scale and is only damped on galactic and 
smaller scales by the nonlinear K-mouflage mechanism. 
See Sec. EH for a discussion and comparison with other 
modified-gravity theories. 


where SM = 47 t<5 (< r)pr 3 / 3. Then, Eq. dSSl) gives for the 
evolution of the normalized radius y : or density contrast 
${< r) = y~ 3 - 1, 


A 2 y 

d(lna) 2 + 


1 dH \ 
~ J+ H 2 ~dt ) 


Ay 

d In a 


+%(i+ei ){y-^ 


1 ) 2 /= 0 . 


(60) 

Again, as in the A-CDM cosmology but in contrast with 
f(R), Dilaton and Symmetron models, the spherical col¬ 
lapse is scale invariant so that the dynamics of differ¬ 
ent mass shells decouple. This applies to the unscreened 
regime, from clusters of galaxies up to the Hubble radius. 


E. Halo mass function 


As usual, we can write the halo mass function 
n(M)AM/M as 


= —t{v) —, 
v ’ M M w v 


with v = 


r{MY 


(61) 


where we used the fact that the linear growing modes 
are scale independent [so that Sl/(t(M) = 6 u/(Ti(M), 
where the subscript “i” refers to the high redshift Z{ where 
the Gaussian initial conditions are defined, far before the 
dark-energy era]. Here a(M) is the root mean square of 
























the linear density contrast at scale M and Sl is the lin¬ 
ear density contrast associated with the nonlinear density 
threshold A m that defines the virialized halos. The map¬ 
ping Sl —> 5 is obtained by solving the spherical collapse 
dynamics m, with the initial condition ?/; = 1 — <5ii/3 at 
a very high redshift z\. Inverting this relation gives the 
linear density threshold Sl that is associated with a given 
nonlinear density threshold S = A m , where the subscript 
“m” denotes that S = y~ 3 —l is the density contrast with 
respect to the mean density of the Universe. 

The scaling variable v directly measures the probabil¬ 
ity of density fluctuations in the Gaussian initial condi¬ 
tions. Then, we take for the scaling function f(v) the fit 
to A-CDM simulations obtained in jJIj], which obeys the 
exponential tail f(y) ~ e~ v / 2 at large v. This means 
that the mass function (IbTI) shows the correct large-mass 
tail, which is set by the Gaussian initial conditions. 


F. Planck masses 


It is interesting to note that, depending on the physical 
process that is considered, one can define several effec¬ 
tive Planck masses. In other words, if we assume Gen¬ 
eral Relativity and measure the reduced Planck mass or 
Newton’s constant from different sets of observations, we 
would obtain different values. This could be used as a 
signature of the modified-gravity theory. 

From Ea. tEBl) . the effective Planck mass that would be 
read from the Friedmann equation, at the background 
level, is 

M 2 (t) - ^Pii 1 ~ e cqo\ 

7 W Pl(Friedmann)UJ ~ - A 2 (t) -’ ^ 

On the other hand, with respect to large-scale den¬ 
sity fluctuations in the cosmological unscreened regime, 
where the Klein-Gordon equation can be linearized over 
the scalar field, the effective Planck mass that would be 
read from the modified Poisson equation (l45j) is 


^(unscreened )(*) ~ A 2 (t)(l + {63) 

On small astrophysical scales, within the highly nonlinear 
screened regime, the effective Planck mass is instead the 
one defined in Eq.flST]), 


M 2 

-^Pl(screened) W = ' (®4) 

The difference between these various definitions is a sig¬ 
nature of the modification of gravity associated with the 
scalar-field models m, as seen from the Jordan frame. 


G. Symmetry of large-scale gravitational clustering 

On large scales, where we neglect shell crossing and 
pressure, the dynamics of gravitational clustering is given 


by Ecis. dSOl) . (l48l) (with p = 0), and (l45|) . It is convenient 
to define the rescaled velocity and metric potential by 


da 


v=— /u, $= — / 


d t 


da 


d t' 


(65) 


where we introduced the linear growth rate 

d In D + 
din a ’ 


( 66 ) 


and to change the time variable from t to In D + . Then, 
the equations of motion write as 


05 

dhiD+ 


+ V ■ [(1 + <S)u] 


<9u 

d In D + 


(u ■ V)u + (k — l)u 


V 2 0 


0, (67) 

-V0, (68) 

k5 , (69) 


where we introduced the time-dependent factor n(t), de¬ 
fined by 


n(t) 


4nSp(l + ei) 

7 d In D + \ 2 

\ dt ) 


30 m 

2/2 


(1 + G). 


(70) 


Therefore, after the change of time coordinate t In D + 
and the rescaling (lC>5l) . the only dependence on cosmol¬ 
ogy that is left in large-scale gravitational clustering is 
encapsulated in the function n(t). This remains valid be¬ 
yond shell crossing but it breaks down on small scales 
where baryonic effects become important and introduce 
new characteristic scales, which cannot be absorbed by 
the change of variables (ESI). Nevertheless, on large scales 
where gravity is the dominant process, this symmetry 
means that all cosmological scenarios with the same func¬ 
tion k(D + ) show the same density and velocity fields 
{<5, u}. In particular, this means that quintessence mod¬ 
els, where only the background dynamics is modified (i.e., 
the Hubble expansion rate H{z)), and modified-gravity 
models or dark-energy models (with dark-energy fluctu¬ 
ations) that only give rise to a modification of Poisson 
equation by a time-dependent Newton’s constant, are 
equivalent with respect to gravitational clustering, if they 
show the same function n(D + ). In the context of A-CDM 
cosmology, this property has also been used to derive ap¬ 
proximate consistency relations satisfied by the matter 
correlation functions [valid at the nonlinear level within 
the approximation where the dependence on cosmology 
of k(D + ) can be neglected] [2jl - f25l . 

In the case of the K-mouflage scenarios, this symmetry 
only holds on large scales (down to cluster scales) where 
the Klein-Gordon equation can be linearized over ip, as 
in Ea. (H51) . On smaller scales (galactic scales and be¬ 
low), higher-order terms over ip become important and 
the nonlinear K-mouflage screening mechanism comes 
into play. Then, the modified Poisson equation no longer 
takes the linear form fi5l) and the symmetry (1701) breaks 
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down. On even smaller scales we actually recover General 
Relativity as $ ~ T N . because the fifth force is screened. 
In hierarchical scenarios, where smaller scales collapse 
first, larger scales are not very sensitive to the details 
of the clustering on smaller scales, while small collapsed 
scales are sensitive to the clustering up to the largest 
scale that has turned nonlinear. Therefore, we expect 
the symmetry GS) to apply to scales that are greater 
than the transition to the K-mouflage screening regime, 
and not to smaller scales (even though we recover Gen¬ 
eral Relativity far inside the nonlinear screening regime). 
Besides, on such small scales non-gravitational baryonic 
effects come into play (such as AGN feedback) and the 
symmetry no longer holds. 

For our purposes in this paper, the formulation (167|) ~ 
fl75D explicitly shows that, from cluster scales up to the 
Hubble scale, K-mouflage cosmologies belong to the same 
family as the A-CDM and quintessence scenarios, with 
respect to matter clustering. The equations that gov¬ 
ern the gravitational dynamics in these rescaled variables 
take the same form, except for a time-dependent function 
k(D + ). However, the shape of this function is similar 
for realistic scenarios (we shall see in Fig. [5] that ei is 
about 2%). Therefore, we can expect that gravitational 
clustering shows the same qualitative properties in these 
cosmologies and only small quantitative deviations. In 
particular, semianalytical methods should work equally 
well for all these cosmologies, and phenomenological ob¬ 
servations, such as the fact that virialized halos are well 
described by Navarro-Frenk-White (NFW) profiles [26] in 
A-CDM cosmology, should remain valid in other cases. 
This justifies our modelization of clusters, described in 
Sec. |V] below, where we treat A-CDM and K-mouflage 
cosmologies in the same manner. 

IV. NUMERICAL RESULTS FOR 
LARGE-SCALE STRUCTURES 

In this paper we consider two simple models for K (x). 
The first one, which we call the “arctan model” in the 
following, is defined by: 

K arc tan(x) = -1 + X + K *[X ~ X* arctan(x/x*)], (71) 
with the low-x expansion 

X-> 0 : A arctan (x) = —1+xH— 7T~‘ 7—(72) 

3 X* 5 x* 

and the choice of parameters 

A* = 10 3 , x* = 10 2 - (73) 

This gives a K-mouflage model that is consistent with 
both Solar System and cosmological constraints (with 
f} = 0.1). For comparison we also consider the model 
used in [ 3 , 3 , which we call the “cubic model” in the 
following, in which K{\) is a low-order polynomial: 

Acubic(x)=-l+X + A'oX m , (74) 


and the choice of the parameters 

m = 3, K () = 1. (75) 

Here Eq. (1711) should not be understood as a perturbative 
expansion around x = 0. It is rather a simple model that 
interpolates between the low-x behavior © and a large-x 
power-law behavior oc x m . The cubic model is consistent 
with cosmological constraints but its form at large neg¬ 
ative Xi — X ^ 1; is n °t consistent with Solar System 
constraints. Therefore, this is an effective model that 
applies to the semiaxis \ > while the large-negative 
domain is left unspecified. This is sufficient for our pur¬ 
poses, since the cosmological background and large-scale 
perturbations correspond to x > 0 and clusters corre¬ 
spond to x > — 1 (more precisely, |x| 1). 

For both models we choose an exponential form for the 
coupling function, 

A(tp) = eW* ipi t with 0 = 0.1 (76) 

We also consider a reference A-CDM model for compar¬ 
ison. 

All the cosmological scenarios are normalized to 
the same background cosmological parameters today, 
{D m o, f2( r )o, Ddeo, Hq}. In addition, we normalize the 
Planck mass (E© to the same value Mp 10 today, as mea¬ 
sured by Solar System experiments. This means that we 
renormalize the Einstein-frame Planck mass by a factor 
Aq, where A 0 = A(z = 0 ), 

A/|i(f) = Mp 10 ’ w h ence Mpi = Mp W A^. (77) 

On the other hand, the matter density power spectrum 
P(k) is normalized to the same value at high redshift, 
when dark energy is subdominant and both cosmologies 
almost coincide. However, these different scenarios do 
not exactly converge in terms of the background expan¬ 
sion rate at a given matter density, because of the differ¬ 
ent high-redshift reduced Planck masses. Therefore, the 
normalization to the same power spectrum for the matter 
density contrast at high z is somewhat arbitrary, since the 
K-mouflage and A-CDM models do not coincide. Never¬ 
theless, for our purposes this is a convenient choice as it 
illustrates how the difference in the gravitational cluster¬ 
ing dynamics that appear at low z, because of the fifth 
force mediated by the scalar field, affect the late-time 
density field. [For other normalization choices, it would 
be difficult to distinguish the effects due to the different 
normalizations at high z, before the dark-energy era and 
when the fifth force was negligible, and to the late-time 
dynamics characterized by different growth rates.] This 
normalization also corresponds to the case where we can 
measure the density contrast field, i.e. the patterns of 
large-scale structures (e.g., the scale associated with the 
nonlinear transition a 2 = 1 ) at high z, independently 
of accurate measures of the background density and ex¬ 
pansion rate. In practice, this normalization ambiguity 
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z 


FIG. 1: Evolution with redshift of the matter and dark-energy 
cosmological parameters fl m (z) and fide (a). We display the 
two K-mouflage models of Ea. dTIII (arctan model, red crosses) 
and Eg. (1711) (cubic model, blue squares), and the reference 
A-CDM universe (black dashed lines). The two scalar-field 
models almost coincide in this figure. 



z 


FIG. 2: Relative deviation of the Hubble expansion rate with 
respect to the A-CDM reference, A H/H = H / Ha-cdm — 1, 
for the same K-mouflage models as in Fig. [T] 


does not appear because one compares each cosmologi¬ 
cal scenario with the data, rather than comparing with 
a theoretical reference cosmology (in particular, the best 
fits associated with different theories will typically have 
slightly different cosmological parameters and expansion 
rates at both z = 0 and z —> oo). 

In the following, we present our results for the choice 
of cosmological parameters today given by D m n = 0.25, 
fl de0 = 0.75, h = 0.70 and cr 8 = 0.7. 



z 


FIG. 3: Effective equation of state parameters Wde(z) (solid 
lines with a divergence and change of sign at z — 3) and 
w v {z ) (dashed lines with a smooth behavior), for the same 
K-mouflage models as in Fig. [T] 


A. Background dynamics 

We consider the density parameters fl m (z) and flde(z) 
in Fig. |T| for the two K-mouflage models defined in equa¬ 
tions CD and (ED and for the reference A-CDM uni¬ 
verse. Since we normalize the density and dark-energy 
parameters to be equal to the ones observed today, all 
models coincide at z = 0 in terms of background quan¬ 
tities. The deviations from the A-CDM scenario are 
slightly greater for the arctan model (I7p than for the 
cubic model d74l) , in agreement with [12}. This is due to 
the fact that K'(x) is slightly smaller in the former case 
over the range of redshifts of interest, z < 6, and that 
deviations from the A-CDM scenario typically scale as 
fi 2 /K', see for instance Ea. (l46l) and Ref.p}. 

As for studies performed in the Einstein frame 3] 
(where the cosmological parameters are normalized by 
their Einstein-frame values today), we find that the dark- 
energy density becomes negative (and subdominant) at 
high redshift. This gives Cl m > 1 at high z for the two 
K-nrouflage models CD and ED (but as i n the A-CDM 
case fl m —^ 1 at high redshift). 

In Fig. [D we consider the relative deviation, 
H(z)/H a-cbm(z) — 1, of the Hubble rate with respect to 
the A-CDM reference. The deviation is slightly larger for 
the cubic model <1731) at 2 ~ 6, but this is only a transient 
effect because at z > 12 the deviation is slightly greater 
for the arctan model ED- as expected. At moderate 
redshifts the Hubble expansion rates differ by less than 
2% between the three cosmologies that we consider here. 
This amplitude is mostly set by our choice of coupling 
constant j3 = 0.1, because as recalled above deviations 
from the A-CDM scenario scale as p 2 /K' and at low z 
we have K' ~ 1. Therefore /3 = 0.1 typically leads to 
percent deviations from the A-CDM scenario. This value 
of (3 (or lower values) is required to satisfy observational 
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constraints on cosmological and Solar System scales [jjj , 
in particular from the expansion rate at the time of Big 
Bang Nucleosynthesis and from the bounds on the cur¬ 
rent time derivative of the gravitational coupling Q pro¬ 
vided by the Lunar Ranging experiment. 

The deviation from the A-CDM reference does not van¬ 
ish at high redshift because the reduced Planck masses 
are different, see Eq. (El). Indeed, in the K-mouflage 
models Mpj(t) becomes time dependent and we choose 
to normalize all scenarios by their Planck mass today 
(when Solar System measurements and laboratory exper¬ 
iments are performed). Note that in studies performed in 
the Einstein frame Q, where the reduced Planck mass is 
constant, one can recover the A-CDM expansion rate at 
both 2 = 0 and at high redshift. However, this requires 
normalizing the matter density today by the conserved 
density p of Ea. dAldl) instead of the Einstein-frame den¬ 
sity p. 

In Fig. [3]we display the effective equation of state pa¬ 
rameter for the dark energy, Wde = Pde/Pde , evaluated 
using Eos. (l29l) and (Iddl) for pde. and pde- For both mod¬ 
els Wde —> — 1 at late times, mimicking the presence of 
a cosmological constant. As in studies performed in the 
Einstein frame, the effective equation of state parameter 
is beyond —1 at low z and changes sign at a moderate 
redshift while going through ±oo (this does not corre¬ 
spond to a singularity in terms of the Hubble rate or 
dark-energy density but to the vanishing and change of 
sign of pde) • We also show in Fig. [3] the equation of state 
parameter w v defined as 


= Pv_ = K 
P<e 2 IK'-K' 


(78) 


where we used Eq. m- In contrast to wde, remains 

negative over 2 < 6 and shows no divergence. The dif¬ 
ference between the behaviors of Wde and w v shows the 
impact of the coupling between the matter and scalar- 
field components. This makes the dark-energy density 
and pressure significantly different from the bare scalar- 
field ones, see Eas. (l29l) and (Iddl) . and can even make pde 
and p v have different signs. 


B. Background scalar field 

We show in Fig. [4] the background values (p and \ 
of the scalar field and of its kinetic term. The scalar 
field (p is negative and its amplitude grows with red¬ 
shift (we chose the normalization (p = 0 at high redshift, 
2 — > oo). We can see that \(3<p/Mpi\ <C 1 until z = 0. 
Thus, the coupling function A(tp) is dominated by its 
low-order terms in the expansion (|8|) and choosing for in¬ 
stance A(p>) = (1 + Pip/nMpi)” , with n > 0, would give 
similar results to the exponential choice (1761) . 

The kinetic term \ decreases with time. It goes to in¬ 
finity at high redshift, 2 —> oo, and we can check that 
at low 2 it is significantly smaller than unity. Then, 



z 



z 


FIG. 4: Upper panel: background scalar field (p as a function 
of redshift. Lower panel: background kinetic term \ as a 
function of redshift. 


the kinetic function K (y) is dominated by its low-order 
terms in the expansion ©• This explains why the two 
K-mouflage models converge at low 2 in the lower panel 
of Fig. m 


C. Modified gravitational potentials, gravitational 
slip and effective Newton’s constant 

As seen in Sec. m both for the background quantities 
and the large-scale perturbative structures, the devia¬ 
tions from the A-CDM universe can be measured by the 
two coefficients e\(t) and e 2 (i), defined in Eas. (H7>l) and 
(1251) . In particular, from Eqs. m and (l44l) the two grav¬ 
itational potentials $ and ’P of the Jordan-frame metric 
m read for large-scale unscreened structures as 

$ = (l + ei)^N, T = (1 — ei)W N , (79) 
and the normalized gravitational slip writes as 
\P - $ 


(80) 
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FIG. 5: Coefficients ei and £ 2 , defined in Egs. (l46l) and (1251) FIG. 6 : Relative drift with redshift of the effective Newton 
for the K-mouflage models, as functions of redshift. constant for the K-mouflage models. 


We show both coefficients ei(f) and e 2 (f) in Fig. [5] We 
can see that they are of order 2% at low z. Again, as can 
be seen from Eq. (l46l) , this amplitude is set by our choice 
/3 = 0.1 (to satisfy observational constraints) as K' ~ 1 
at low z and deviations from the A-CDM reference then 
scale as f3 2 . This also sets the amplitude of the gravita¬ 
tional slip 77 . At high 2 the coefficients et(t) go to zero, as 
K' goes to infinity and we enter a cosmological nonlinear 
screening regime that also ensures that the dark-energy 
component becomes subdominant at early times. This 
decrease of appears faster for the cubic model, be¬ 
cause of its stronger growth of K'(x) at large positive %. 
As noticed above, this means that departures from the 
A-CDM scenario are greater for the arctan model than 
for the cubic model (with our choice of parameters). 

The Jordan-frame coefficient ei is always positive and 
the gravitational slip 77 defined in Eg. ((501) is negative. 
From Eg. (l25l) the coefficient e 2 also reads as 



z 


FIG. 7: Relative deviation, D+/D+ a-cdm — 1, of the linear 
growing mode D+ from the A-CDM reference. 


e . ( t ) = P® dip = P® (81) 

21 ‘ Mpi d In a M Pl H(t)dt’ [ 

which is negative from Ea. (155l) and of order /3 2 /K'. 

In Fig. [G] we show the evolution with redshift of the 
effective Newton constant, defined from Eq. (El) as 


Git) = Go 


A 2 (t) 

Al 


(82) 


Because of the dependence of the effective Newton cou¬ 
pling strength on the background value of the scalar field, 
G is a few percent higher at z ~ 6 than today. 


D. Linear theory 



z 


In Fig. □ we show the relative deviation, 
D + /ZI + a-cdm — 1 , of the linear growing mode, 
obtained by solving equation (l55l) . from the A-CDM 


FIG. 8: Relative deviation, f / f acdm — 1, of the linear growth 
rate, f(z) = dlnD+/dlna, from the A-CDM reference. 
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reference case, and the linear growth rates f{z). Again, 
the relative deviation of the growing mode is greater for 
the arctan model CD than for the cubic model m , 
because of the lower value of K' over relevant redshifts, 
see Ea. (l46l) for the coefficient that modifies the linear 
growing mode equation (15511 . All linear growing modes 
converge at high redshift, despite the slightly different 
Planck masses and Hubble expansion rates. Indeed, 
far in the early matter-dominated era we recover an 
Einstein-de Sitter cosmology and the Hubble term in the 
parenthesis in Ea. (l55l) converges to H~ 2 dH/dt —► —3/2. 
Moreover, the factor goes to zero because of the 
nonlinear K-mouflage screening mechanism, see Eq. (1461) , 
as at high redshift \ and K' become large. This large-AT' 
behavior is also required to ensure that the background 
dark-energy density becomes subdominant. 


We show the relative deviation of the linear growth 
rates f(z) in Fig. [8] Overall, f(z) is greater for the K- 
mouflage scenarios, in agreement with the higher value of 
the linear growing mode shown in Fig. [7] The deviation 
is again of the order of a few percent. The deviation for 
the linear modes D + shown in Fig. [7] could reach 10% 
at z = 0 because of the cumulative effect due to the 
integration over time. The growth rates f(z) converge to 
unity at high redshift but we can see that there remains 
a noticeable difference between the K-mouflage models 
and the A-CDM reference up to z ~ 6. 


This rather slow decrease of the deviations from the A- 
CDM reference at higher redshift is a characteristic sig¬ 
nature of K-mouflage models, as many other modified- 
gravity models, such as f(R) theories, lead to a faster 
convergence to the A-CDM scenario at z > 2. This is 
related to the fact that in the linear sub-horizon regime 
the K-mouflage effects are scale independent, as the fac¬ 
tors ei(t) and £ 2 (t) only depend on time. In contrast, in 
f(R) theories or Dilaton models, the factor e(k,t) that 
appears in the modified Euler or Poisson equations, or in 
the evolution equation for the linear density modes, takes 
the form e(k, t) oc /3 2 k 2 /(a 2 m 2 + k 2 ) 1 with a characteristic 
physical scale 27r /m beyond which the theory converges 
to General Relativity. At high redshift this scale typ¬ 
ically goes to zero, so that at a fixed physical (or also 
comoving) scale, deviations from the A-CDM scenario 
vanish because the coupling /3 decreases or one enters 
the unmodified regime beyond 27r/m. In the K-mouflage 
models that we consider in this paper, because there is 
no such characteristic scale the convergence to General 
Relativity is only due to the vanishing of the effective 
coupling strength f3 2 /K', with /3 being constant (in our 
case) and K' increasing at high z because of the nonlin¬ 
ear K-mouflage mechanism itself. However, this decrease 
of f3 2 /K' at high z is rather slow for generic kinetic func¬ 
tions K (y), as seen from the curve obtained for ei(t) in 
Fig. [3] for the simple cubic model. 



z 


FIG. 9: Linear density contrast threshold ^n(A) associated 
with the nonlinear density contrast A c = 200 with respect to 
the critical density p C rit, for the K-mouflage models and the 
A-CDM reference. 
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FIG. 10: Relative deviation, n(M) /riA-CDM(M) — 1, of the 
halo mass function of the K-mouflage models from the A- 
CDM reference, at 2 = 0 (solid lines) and 2 = 2 (dashed 
lines). Halos are defined by the density contrast A c = 200 
with respect to the critical density. 


E. Halo mass function 

Solving the spherical collapse equation (l60l) . we can 
compute the linear density contrast threshold Sl(z) that 
corresponds to a nonlinear density contrast of = 

A m , where A m is the nonlinear threshold that we choose 
to define halos. As discussed in [g), we are not interested 
in Sl(z) at the observation redshift, because it is not an 
observable quantity. Instead, we wish to evaluate the 
linear threshold Sl { , at a given high redshift z\, that is 
required to produce at later time z a nonlinear density 
contrast A m . In other words, we want to estimate the 
initial density fluctuation associated with a given nonlin¬ 
ear density contrast at the observed redshift. To avoid 
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the introduction of an arbitrary high redshift Zi, follow¬ 
ing what it is done in [s| and usual practice, we trans¬ 
late all the initial thresholds Sl, to redshift 2 by multi¬ 
plying them by H+a-cdm^)/£>+A-CDM( zi) [instead of 
using D + (z)/D + (zj), i.e. the linear growing mode as¬ 
sociated with each K-mouflage scenario], and we denote 
this quantity the linear density contrast threshold 
This common multiplicative factor enables a meaningful 
comparison between the different scenarios. 

In Fig [9] we display <5l(a) when we define halos by a 
constant density contrast threshold A c = 200 with re¬ 
spect to the critical density p cr it- This corresponds to a 
density contrast with respect to the mean density p of 
A m = A c /f2 m (.z). We choose a constant A c rather than 
A m because observational cluster surveys usually define 
cluster halos by a constant overdensity with respect to 
the critical density p cr it . At high redshift both definitions 
coincide, as fl m —> 1 , while at low redshift or far in the 
dark-energy dominated era one can argue that A c makes 
more physical sense. Indeed, during an exponential accel¬ 
eration phase the growth of large-scale structures freezes 
out and one obtains isolated halos among increasingly 
large voids. Then, the mean universe density p decreases 
as a -3 , following the dilution due to the expansion, and 
it does not correspond to the typical density of halos (nor 
voids). In contrast, we can assume that the isolated halos 
keep a roughly constant physical radius and density, like 
the critical density p cr it (in an exponential phase where 
the Hubble rate is constant), so that it is more meaning¬ 
ful to express halo densities in terms of /0 C rit- 

Both K-mouflage models accelerate the growth of 
large-scale structures as compared with the A-CDM ref¬ 
erence, as seen from the linear growing modes of Fig. [7] 
Therefore, we find in Fig]9] that a smaller linear density 
contrast is required to reach the same nonlinear overden¬ 
sity of 200 (with respect to the critical density). Again, 
the deviation from the A-CDM prediction is greater for 
the arctan model CD. 

From the linear threshold displayed in Figd we ob¬ 
tain the halo mass function as in Ea. lRJlT) (note that 
v = S L /cr = 5 L (a)/ct(a) = SLi/cri). In Fig [10] we show 
the relative deviation n(M) /tia-cdm(M) — 1 of the halo 
mass function from the A-CDM reference case, at z = 0 
and z = 2. As for the case of S L ( a), since the scalar 
field enhances gravitational clustering we find that the 
mass functions for the two K-nrouflage models are higher 
than the A-CDM reference in the high-mass tail, with 
the greater deviation obtained for the arctan model CD- 
As usual, the deviation increases at high mass because 
the exponential falloff amplifies the sensitivity to slight 
departures of the growth of structures. (The deviation 
becomes slightly negative at low mass because all mass 
functions obey the same normalization to unity: there 
cannot be more matter in halos than the matter content 
of the universe.) 

At fixed mass, M ~ 10 14 /i _ 1 M 0 , the deviation from 
the A-CDM reference is greater at z = 2 than at z = 0, 
despite the difference in linear density thresholds being 
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FIG. 11: Upper panel: relative deviation of the matter density 
power spectrum from the A-CDM reference, at 2 = 0 (solid 
lines) and 2 = 2 (dashed lines), for the K-mouflage models. 
Lower panel: relative deviation of the matter density correla¬ 
tion function from the A-CDM reference. 

lower, as seen in Fig. [9] This is because at fixed mass we 
are further into the rare high-mass tail, which amplifies 
the dependence on the linear density threshold and more 
than compensates the slow convergence between the K- 
mouflage and A-CDM thresholds. 


F. Matter density power spectrum and correlation 
function 

We show in Fig. [TT] the matter density power spec¬ 
tra and correlation functions at z = 0 and 2 = 2. The 
computation of this power spectrum P(k) combines per¬ 
turbation theory up to one-loop order with a halo model, 
as described in Q in the case of Einstein-frame stud¬ 
ies and following the approach introduced in [27] . The 
two-point correlation function £(x) is obtained from the 
Fourier transform of the power spectrum. On large scales 
we obtain the same constant relative deviation for the 
power spectra and the correlation functions, as the lin¬ 
ear growing modes D + (t) are scale independent in both 
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K-mouflage and A-CDM cosmologies (much below the 
horizon). This deviation of 20% is consistent with the 
deviation of 10% obtained for the linear growing modes 
in Fig. 0 The deviation from the A-CDM reference grows 
on mildly nonlinear scales, as nonlinearities amplify the 
effects of the fifth force. This is related to the increase 
of the large-mass tail of the halo mass function shown in 
Fig. [TUI because on these scales the power spectrum and 
the correlation function probe the formation of massive 
halos, as can be clearly seen in a halo model approach. 
The deviation decreases on smaller scales because the 
power spectrum now probes the inner regions of halos and 
we assume similar NFW profiles [26} for all cosmologies 
(but this regime shows a greater theoretical inaccuracy 
and numerical simulations would be need to measure the 
impact of the modified gravity on small highly nonlinear 
scales and halo profiles). However, in the nonlinear range 
shown in Fig.[lTJ the impact of changes to the mass func¬ 
tion is greater than that of halo profiles, see also Ref. [28jj| 
for a detailed study of these various contributions. 

The deviation from the A-CDM reference slowly de¬ 
creases at high z, as the fifth force mediated by the scalar 
field becomes negligible (as seen from the vanishing of the 
key factor /3 2 /K' as K' —> oo). This decrease of the devi¬ 
ations of large-scale clustering from the A-CDM reference 
is slower than what is found in many other modified- 
gravity scenarios, such as the f(R) theories, and is char¬ 
acteristic of these K-mouflage models. 


V. CLUSTERS OF GALAXIES 

To go beyond background quantities and the large- 
scale perturbative regime, we investigate in this section 
the impact of K-mouflage scenarios on the largest col¬ 
lapsed structures that we observe, that is, clusters of 
galaxies. This provides another probe of modified-gravity 
models, which is complementary with background and 
perturbative studies, as it corresponds to the nonlinear 
regime of the matter density field and to the well-defined 
objects measured in actual surveys. 

For our purposes, clusters present two advantages as 
compared with galaxies. First, they are unscreened ob¬ 
jects [If, so that the impact of the modification to gravity 
is very simple and corresponds to a time-dependent ef¬ 
fective Newton constant (we shall check that this holds 
down to the cluster cores in Sec. IVBl below). There¬ 
fore, one does not expect dramatic qualitative changes 
from the A-CDM case, and the same semi-quantitative 
models can be applied to both K-mouflage and A-CDM 
cosmologies. This is also illustrated by the symmetry 
described in Sec. IIII Gl which shows that in this un¬ 
screened regime, from the point of view of nonlinear grav¬ 
itational clustering, the A-CDM cosmology, quintessence 
models, and K-mouflage scenarios, belong to the same 
class. They obey the same equations of motion (lUTll - flUSil . 
with only slightly different time-dependent functions n(t) 
from Eq. HDD. Second, at first order clusters can be de- 
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FIG. 12: Mass-concentration relation for NFW halos, for the 
K-mouflage models and the A-CDM reference, at 2 = 0.37. 
The black points (with their error bars) are observational 
measures taken from [3Cj]. 


scribed by simple physical laws, such as hydrostatic equi¬ 
librium for the gas profile and Bremsstrahlung emission 
for the X-ray luminosity, giving rise to the so-called clus¬ 
ter “scaling laws” [29|. This is especially true for the 
most massive clusters that we focus on. 

In contrast, galaxies probe the transition from the un¬ 
screened to screened regimes and also involve many com¬ 
plex astrophysical phenomena, such as cooling processes, 
star formation, supernovae and AGN winds and feed¬ 
back. Therefore, although they would be very interest¬ 
ing probes they are much more difficult to model, both 
for the modified-gravity sector and for the usual galaxy 
formation processes that also appear in the A-CDM cos¬ 
mology. 

In this paper, our goal is to estimate the magnitude of 
the impact of K-mouflage models on clusters of galaxies 
rather than building a very accurate description of clus¬ 
ters. Therefore, we consider the simplest possible mod¬ 
elling of clusters with only few physical parameters. This 
may not provide the highest-accuracy cluster model, but 
we can expect that it captures the main physical pro¬ 
cesses and provides a robust estimate of the impact of 
modifications to gravity. Moreover, we check that our 
predictions show a reasonable agreement with observa¬ 
tions. 


A. Halo density profiles 


To study the effects of K-mouflage scenarios on clusters 
of galaxies we need their dark matter profile. Because in 
the unscreened regime gravitational clustering proceeds 
in the same fashion in K-mouflage and A-CDM cosmolo¬ 
gies, as described in Sec. IIII Gl we assume in all cases 
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NFW profiles [26| . 


PdmM 


_ Ps _ 

(r/r s )(l +r/r s ) 2 ' 


(83) 


This profile is characterized by a scaling radius and den¬ 
sity, respectively r s and p Sl which can be expressed in 
terms of the concentration parameter c = R/± c /r s . Here 
i?A c is the radius such that the mean density within i?A c 
is A c times the critical density, p{< Ra c ) = A c p cr it, as 
we again define the extension of the halos by an overden¬ 
sity threshold with respect to the critical density. From 
the definition of c, is possible to express p s as 

A . c 3 

Ps = Pcrit “X - 1 JZ ; \ J7~, ' (84) 

3 in(l + c) — c/(l + c) 


which can be inverted to give c as a function of p s . 

To consider the effects of the presence of the scalar 
field on the shape of the dark matter profile we consider 
a simple model for the halo concentration. We assume 
that halos of mass M typically form at a redsliift zt(M) 
determined by 


<r(q,z { ) = a u (85) 

where q = (SM/Anpo) 1 ^ 3 is the halo Lagrangian radius 
and <7f is a free parameter, and that the density of the 
newly-formed halo is proportional to p c r it(zf), 

Ps(M) = Af Pcrit(Zf), (86) 

with Af a second free parameter. Equation (1851) means 
that halos of a given mass typically form when density 
fluctuations at this mass scale reach the nonlinear regime, 
while equation (1551) assumes that the core of the clus¬ 
ter keeps a roughly constant density after its formation, 
which is set by the critical density at the formation time. 
As discussed in Sec. IIVE1 we choose to rescale p s in terms 
of the critical density rather than the mean density at 
redshift Zf because the former is more physical at late 
times (whereas they coincide at high redshift) and it also 
corresponds to our definition of halos. Next, using equa¬ 
tion (1841) . we compute c(M) and we define the dark mat¬ 
ter density profile using (1551) . 

In Fig. [T5] we display the mass-concentration relation 
obtained with the choice of parameters at = 0.2 and 
Af = 500 (halos are again defined by A c = 200). As 
is well known P, the concentration c (and the scaling 
density p s ) is larger for smaller mass, because in hierar¬ 
chical scenarios smaller mass scales turned nonlinear at 
higher redshift when the critical (and the mean) density 
of the Universe was greater. This is of course consistent 
with our model (1551) - (1551) . We compare these results with 
the mass-concentration relation obtained by [30], from 
the analysis of 19 X-ray selected galaxy clusters from 
the Cluster Lensing and Supernova Survey with Hubble 
(CLASH), with a mean redshift z ~ 0.37. We can see 
that reasonable choices of the parameters erf and Af (we 
naturally expect at < 1 and Af > 200) allow us to obtain 


a reasonable match to observations. This suggests that 
this simple modelling captures the main features of the 
gravitational formation of halos. Therefore, we do not 
consider here more sophisticated models, which involve 
the past accretion history and merging trees of virialized 
halos. These could provide more accurate modelling, at 
the price of additional complexity (and often additional 
parameters), but it is not clear if their estimate of the 
dependence on the underlying gravity theory would be 
much more accurate. Such studies are left for future 
works, where N-body simulations would be needed to 
check detailed models. 

As expected, we find a small increase of the concentra¬ 
tion c(M) in the K-mouflage models, as compared with 
the A-CDM reference. This is due to the faster growth of 
gravitational clustering, which implies a slightly greater 
scaling density p s (M). However, we can see that the 
effect is rather modest. 


B. Clusters are not screened 

As noticed in ®], clusters are unscreened, and the 
Klein-Gordon equation (1551) can be kept at the linear level 
over the fluctuations of the scalar field, as in Eq. (1551) . In 
this section we check that this property extends far inside 
the cluster profile. 

In the small-scale static limit, which corresponds for 
instance to high-density astrophysical objects, the Klein- 
Gordon equation (IA3D becomes 


V r ■ (V r (/7 K') = p, (87) 

(rM pi 

where r is the Jordan-frame physical coordinate and we 
assumed 5p~ p (i.e., S 1). As compared with Eg. (1551) . 
here we do not make the approximation I\' ~ K'. In¬ 
stead we consider the small-scale regime where \ < 

0. For a spherically symmetric halo, using the Stokes 
theorem this gives 


d <P K , = /3 AM (r ) 
dr c 2 Mpi47rr 2 


( 88 ) 


As in [7>, 13, [HJ (but in Jordan-frame coordinates), defin¬ 
ing the “K-mouflage screening radius” Rk by 


Rk{M) = 


f3A 2 M 

AitcMpiM 2 


1/2 


(89) 


where M = M(R) is the total mass of the object of radius 
R, and introducing the rescaled dimensionless variables 
x = t/Rk , m(x) = M(< r)/M , 4>(x) = tp(r)/ipK , with 
Pk = M 2 Rk/cA, the integrated Klein-Gordon equation 
(1551) reads as 



m(x) 

9 
X z 


with 


X = 


1 

2 



(90) 
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FIG. 13: Upper panel: scalar-field radial profile, 5p{r) = 
ip{r) — i p, for halos of mass 10 15 1i _1 Mq (upper solid lines) 
and 10 13 /i _1 Mq (lower dashed lines). The scalar-field fluctu¬ 
ation is negative as the scalar field ip is minimum at the center 
of the halo. Lower panel: radial profile of the “kinetic energy” 
—x(r) °f the scalar field, for halos of mass 10 15 h~ 1 AL@ (upper 
solid lines) and 10 13 /i _1 Mq (lower dashed lines). Here X < 0 
because we consider the static limit, which is dominated by 
spatial gradients. The arctan and cubic K-mouflage models 
give almost identical results in this figure. 

As pointed out in [?J, in the small-scale static regime we 
have x < 0 whereas the cosmological background value 
satisfies x > 0. Using that M 4 ~ pdeO is roughly the 
dark-energy density today and that A ~ 1 , we obtain 


Rk{M) ~ a -3470 a.u. ~ w ——-0.12 h 4 Mpc. 

^ ’ y 1 Mg Y 10 14 M 0 1 

(91) 

The first equality shows that the Solar System is screened 
by the Sun, which allows these K-mouflage scenarios to 
satisfy Solar System constraints 0. On the other hand, 
for /3 = 0.1 the K-mouflage screening radius of a cluster 
of mass 10 14 Mq is Rk — 0.04ft _1 Mpc. This is much 
smaller than the radius of the cluster, which means that 
most of the cluster is unscreened. Moreover, as we move 
inside the halo the enclosed mass M(< r) decreases, 


which further delays the onset of K-mouflage screening. 
When |x| -C 1 we have K' ~ 1 and we obtain 

Ixl < 1 = X(r) - > ( 92 ) 

which gives inside the halo 

|x| <1, r<R : x{r) ~ (ftA 2 A c (< r)^j . 

(93) 

We show the radial profiles of — x(r) in the lower panel in 
Fig. H3] for M = 10 l 5 h- 4 M Q and 10 ls h~ 4 M Q . In both 
cases we can check that |x| -C 1 over the full halo profile 
[as seen from Eci. flO . \{r) goes to a finite value at r —> 0 
for NFW density profiles, because p oc 1/r and m oc x 2 in 
the core]. This means that clusters are not screened and 
that we can use the background value K' for the kinetic 
function. In fact, at low z we also have x <C 1 and 
K' ~ 1, so that the kinetic function is dominated by the 
low-order terms in the expansion 0 and the results are 
not very sensitive to the precise nonlinear form of K (x). 
Then, the Klein-Gordon equation can be linearized in the 
scalar field as in Ea. flTfl) and the gravitational potential 
$ that governs the dynamics of matter is again given by 

Eq.gg). 

From the kinetic factor x(r) we can obtain the radial 
profile of the scalar field, ip(r), by integrating d<j>/dx = 
\J— 2y and using <p = (px<f>. The boundary condition is 
ip —> <p at infinity. We show the radial profile of the fluc¬ 
tuation Sip = ip — ip in the lower panel in Fig. 1131 We can 
check from the comparison with Fig. [4] that \5ip\ <C \ip\, 
in agreement with the scaling Sip/ip ~ (aH/ck) 2 6 p/p ob¬ 
tained from Ea. (IA35l) and tp/Mp\ ~ ft/K'. In particular, 
this explicitly shows that the coupling function A(ip) re¬ 
mains dominated by its low-order terms in the expansion 
®, both for the background and for large-scale struc¬ 
tures such as clusters of galaxies. 

The magnitude of Sip can also be read from the modi¬ 
fied Poisson equations (TTTT) . In realistic models the fifth 
force should not have a magnitude greater than the New¬ 
tonian force, which implies |5A/A| ~ \ftSip/Mp\\ < ||, 
whence \/35ip/Mpi\ < 10 -5 . 

C. Impact of nonlinear substructures 

The Klein-Gordon equation (l87l) that determines the 
scalar field ip is nonlinear, because of the factor K'(x). 
This means that substructures could have a strong im¬ 
pact on the solution ip(r) as there is no longer a linear 
superposition property and the solution obtained for the 
averaged halo profile is not identical to the average of the 
exact solutions obtained by taking into account substruc¬ 
tures. In this section, we check that this nonlinearity 
does not play a significant role and does not invalidate 
our approach described in Sec. IV Bl 

First, we note that for an object that is exactly spher¬ 
ically symmetric, the integrated Klein-Gordon equation 
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(l88l) holds and the scalar-field profile only depends on the 
integrated mass M ( r ) within radius r. This smoothes out 
radial substructures. However, in practice clusters are 
not exactly spherically symmetric, and individual cluster 
galaxies also break any overall spherical symmetry. We 
have seen in Sec. IV Bl that clusters are unscreened as the 
kinetic factor, x c ius(r), associated with the mean cluster 
density profile (1831) . is much smaller than unity. Then, 
if galactic halos only form a small fraction of the total 
cluster volume, throughout most of the cluster volume 
we can linearize Eg. (1571) . as in Eg . (1731) . which gives 


unscreened region: V 2 </? 


PA 

K'c 2 M?i P ' 


(94) 


where K' is the background value, with K' ~ 1 at low z. 
This equation breaks down around each cluster galaxy, 
where the high matter density, which is much greater 
than the NFW mean density pdm(^) of Eg. (1531) at that 
radius, makes the scalar field enter the nonlinear screen¬ 
ing regime. Thus, around each galaxy “i”, i = 1, ..,N ga i, 
we must cut a patch Vxi where the equation (1571) must 
be replaced by the fully nonlinear equation (1571) . By def¬ 
inition, the volume Vk% is given by the K-mouflage ra¬ 
dius Rxi of the galaxy. In practice, if Rxi -C -Rcius, we 
can build an approximate solution by patching the so¬ 
lutions within each galaxy volume Vk% with the global 
solution (1571) that holds in between galaxies. Around 
each galaxy, using an approximate spherical symmetry 
around the galaxy center, we obtain the local profile by 
solving Eq. (1551) , where M is replaced by the galaxy mass 
mga.i(r), and the boundary condition at Rki is approx¬ 
imated as a constant obtained from the global solution 

(El). 

This scenario holds provided the regions Vk, do not 
extend far beyond the galaxy volumes V) (defined for in¬ 
stance by their stellar content or by the region where 
matter is gravitationally bound to the galaxy) and do not 
cover most of the cluster volume. From Eq. ( 1551) we have 
Rk (wigai) oc m gal . Defining the mass function n(m)dm 
of the cluster galaxies, the total volume built by the non¬ 
linear regions Vxi reads as 

du d?) 

VK-gai = / dm— I4r(m) oc / dm—m 3/2 . (95) 
J o dm J o din 


The mass function of the cluster galaxies or of dark mat¬ 
ter subhalos is typically a power law at law mass with 
an exponential cutoff at high mass 32]. In any case, 
the integral f dm(dn/dm)m = M ga \ is necessarily finite 
and equal to the total mass associated with the galaxies, 
which is smaller than the total cluster mass. Therefore, 
the integral (1551) converges at low mass and is dominated 
by the galaxies around the knee of the galaxy multiplicity 
function, which typically corresponds to M ~ 10 12 Mq 
From Eq. (1511) we have Rk( 10 12 M q ) ~ 4 h 1 kpc, with 
/3 = 0.1, which gives a volume fraction of the order of 
(0.004) 3 ~ 5 x 10 -8 . Even if we have ~ 20 such galaxies 
in the cluster, this only makes a fraction of order 10~ 6 of 
the cluster volume. 


Moreover, we can see that Rxi is typically smaller than 
the actual radius Ri of the galaxy (by a factor of a few). 
In the A-CDM cosmology itself, the analysis of the hot 
gas that makes most of the intracluster medium and gives 
rise to the cluster X-ray luminosity (based on hydrostatic 
equilibrium and scaling laws) only applies outside of the 
cluster galaxies, where cooling and star formation pro¬ 
cesses play a major role. Therefore, the nonlinearities 
of the Klein-Gordon equation (1571) do not bring further 
restrictions as compared with the standard A-CDM case, 
as they are “hidden” within the galaxies, and the impact 
of the fifth force on the intracluster medium can be ob¬ 
tained from the linearized equation (1571) associated with 
the unscreened regime. They do not modify global prop¬ 
erties either, such as the cluster correlation function, as 
the dynamics and formation of the clusters remain gov¬ 
erned by the linearized Klein-Gordon equation (1941) . 


D. Hydrostatic equilibrium 


From equation (1751) , the equation of hydrostatic equi¬ 
librium for the gas density p g and pressure p g reads as 


V4> = V 



VjJg 

Pg 


(96) 


where we used Eqs. m and (1441) . This explicitly shows 
how the pressure gradient is amplified, at fixed density 
profile, by the fifth force. Assuming spherical symmetry 
this leads to 


d Pg 
dr 


-PgiX + e l) 


(1T \ 
dr 


(97) 


where we used the fact that clusters are unscreened, so 
that K' ~ K' and the Klein-Gordon equation can again 
be linearized as in Eg. (1571) . In the A-CDM cosmology 
we simply have e\ = 0. To obtain the gas profile from 
Eq. ( 1571) we also need an equation of state that gives the 
pressure as a function of the gas density or temperature. 
We consider an isothermal gas with p g = PgksTg/(pm p ), 
where ks is Boltzmann’s constant and pm p is the mean 
molecular weight of the gas. This yields the gas density 
profile 


p g (r) oc 


(98) 


where the Newtonian potential Tn is fixed by the dark 
matter profile. 

To evaluate T g we assume that the gas temperature is 
proportional to the mean value of the dark matter “tem¬ 
perature”, Tdm(^), which we define from the velocity dis¬ 
persion (Tq(t') as 

feJb M (r) = pm p al, M {r). (99) 


The dark matter being collisionless it is not described by 
a thermodynamical temperature. However, we can ex¬ 
pect the virialization processes associated with the for¬ 
mation of the halo to scale in the same fashion for the 
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FIG. 14: Upper panel: density profiles for a cluster of mass 
M200C = 5 x 10 14 /i -1 Mq. The upper solid lines refer to the 
dark matter density profiles and the lower dotted lines to the 
gas density profiles. The K-mouflage models and the A-CDM 
reference cannot be distinguished in this figure. Lower panel: 
relative deviation from the A-CDM reference of the dark mat¬ 
ter (solid lines) and gas (dotted lines) density profiles. 


dark matter, as measured by its velocity dispersion, and 
for the gas, as measured by its temperature. In particu¬ 
lar, the dark matter velocity dispersion obeys the Jeans 
equation, which can be written at equilibrium as fjjt 

J_d[»OM^M] = _J£ = _ (1+<i) l»N (100) 

Pdm dr dr dr 

For a given dark matter halo profile, set by the NFW 
profile (1551) and concentration parameter c(M), the Jeans 
equation (11001) determines the velocity dispersion profile 
(jpM (r), whence the effective dark matter temperature 
^dm (0 defined in Eg. (1551) . Then, we set the gas temper¬ 
ature T g as 


f, }_r _ J_ /dr4 tt r 2 p DM (r)T DM (r) 

9 p g DM fig f dr47rr 2 p DM (r) 


( 101 ) 


clusters of galaxies (34(. In other words, we assume that 
the kinetic and thermal energies of the dark matter and 
the gas are proportional, because they are generated by 
the same process, the formation and virialization of the 
halo. 

Finally, to fully specify the gas density profile (1551) we 
normalize it as 


r R *c o, 

M g = / dr 4nr 2 p g (r) = —— M D m- (102) 

Jo “DM 

Thus, we consider that the baryon and dark matter mass 
fractions in clusters are given by the cosmological abun¬ 
dance. This assumes that there is no significant redistri¬ 
bution and segregation of matter on scales greater than 
cluster radii, which should be a reasonable approximation 
for massive clusters. 

Therefore, in terms of the intracluster medium, the dif¬ 
ferences between the K-mouflage and A-CDM scenarios 
only arise through three effects, in our framework. First, 
the dark matter profiles (1551) are slightly different because 
of the small change of the halo concentration shown in 
Fig. [12] which comes from the different growth rates of 
large-scale structures. Second, the equation of hydro¬ 
static equilibrium (1551) is modified by the factor e\(t), 
which corresponds to the amplification of gravity by the 
fifth force in the unscreened regime. This implies slightly 
greater pressure gradients for the gas. Third, the gas 
(and dark matter) temperature itself is also amplified by 
the same factor (1 + ei), at a fixed dark matter profile, 
because it also arises from the gravitational collapse, see 
Eas. dlOOD and (11011) . The second and third effects com¬ 
pensate in terms of the gas density profile, as the greater 
potential depth is balanced by the greater gas temper¬ 
ature, so that we can expect rather modest deviations 
between the different cosmological scenarios. 

In Fig. [KjJ we show the dark matter and gas density 
profiles for a cluster of mass M = 5 x 10 14 h~ 1 M Q at z = 
0, for the K-mouflage models and the A-CDM reference. 
The presence of the scalar field makes the density profiles 
more compact, in agreement with Fig. 1121 As expected, 
the deviations from the A-CDM reference are of the order 
of a few percent. 


E. Scaling laws 


From the gas density profile p g (r ) and temperature T g , 
we obtain the X-ray cluster luminosity within radius R 


as 


r R „ 

L x (< R) = 4:7re x {T g ) / n 2 g {r)r 2 dr , (103) 

Jo 


where fi g is a free parameter that we fix to be equal to 0.6, where n g (r) = p g (r)/pm v is the cluster gas density and 
which is of the same order as the values used in studies of ex (T) is the X-ray emissivity, which can be expressed in 




















20 


terms of the temperature as [32j 

4 — 9Y / T \ - 1 / 2 

e x (T) = 4.836 x HT 24 ^-^ {^y) erg,-W 

(104) 

Here Y is the Helium mass fraction, and /i = 2/(4 — Y), 
n e /n g = (2 - Y )/(4 - Y), ( n H + 4 n H e)/n g = 2/(4 - Y), 
where we assume complete ionization. This applies to 
high temperatures of order 1 keV and above, where the 
X-ray emissivity is dominated by Bremsstrahlung. Equa¬ 
tion (11031) . with the emissivity (11041) . gives the total bolo- 
metric luminosity. In practice, one measures the radia¬ 
tion from X-ray clusters within finite frequency bands. 
Therefore, we also define the luminosity within frequency 
bands, denoted for instance as bands “A = [*q 4 , , 

a B = \vf,v 2 S ]”,-, by 

L x a(< R) = L x {< R) U-^t/kaTg _ e ~hv£/k B T g \ 

(1 °5) 

Observational studies often measure the X-ray proper¬ 
ties of galaxy clusters within a radius Rx that is smaller 
than -R2000 because the luminosity scales as the squared 
density [see Eg. (11031) ] so that inner high-density regions 
are easier to measure. A popular choice is the radius set 
by the density threshold A c = 500 with respect to the 
critical density. In the following, keeping our definition 
of halos by the threshold A c = 200 as in Figs. [TUI and ITS1- 
M we use the density profile obtained from Ea. (l98l) and 
displayed in Fig. [ll]to compute X-ray properties within 
Rx defined by A c = 500 (hence Rx < -Rhaio)- 

To avoid the complications due to the internal struc¬ 
tures of the clusters (presence of massive galaxies in 
the center, importance of dissipative processes, cooling 
cores,...) and also to follow the observational proce¬ 
dures, we define a core radius R core outside of which we 
evaluate the quantities of interest. As in many observa¬ 
tional analyses, we simply define R c ore as a fixed frac¬ 
tion of the cluster radius Rx (as defined by the thresh¬ 
old A c = 500 with respect to the critical density), with 
-Score = fcoieRx and typically / core ~ 0.15. Then, we 
obtain for instance the luminosity in the outer cluster 
shells, between R coie < r < Rx, as 

Lxa no—core — Lxa{< Rx) — L X a{< R core ). (106) 

In Figs. [15] and [16] we show, respectively, the 

M 50 O 12 ~ 7500c —nocore and T 50 oc —nocore — Lxa —500c—nocore 

relations compared to observations of clusters of galaxies 
in the X-ray, with the choice of parameter / core = 0.15 
and the frequency “A band” [0.1 - 2.4] keV. For the M—T 
relation we obtain a good agreement with observations 
while our prediction for the slope of the T — L relation 
is too shallow. This is a well-known problem associated 
with a noticeable breakdown of the naive “scaling laws” 
for the X-ray luminosity, especially for small clusters [35|. 
This is usually explained by a decrease of the gas frac¬ 
tion and a greater importance of nonthermal effects, or 
departures from hydrostatic equilibrium, in small clus¬ 
ters. However, because our goal is only to estimate the 


> 
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FIG. 15: Upper panel: mass-temperature relation for the K- 
mouflage models and the A-CDM reference, at z = 0.0288 
(lower curves) and z = 0.451 (upper curves). The data points 
are taken from observations made by [3(| (in green), [37| (in 
magenta) and [3(§] (in brown), with clusters in the redshift 
range 0.0288 < z < 0.451. Lower panel: relative deviation 
of the cluster mass-temperature relation from the A-CDM 
reference, at z = 0.0288. 


magnitude of the effects due to modifications of gravity, 
we do not try to build a more accurate and more com¬ 
plex model in this paper. Moreover, our simple model is 
sufficient to recover the typical X-ray luminosity in the 
range 4 < T < 15 keV, which corresponds to massive 
bright clusters. 

At fixed mass, the temperature in the K-mouflage sce¬ 
narios is greater than in the A-CDM reference by about 
2%. This is mostly set by the factor ei, which is about 
2% as seen in Fig. [5] Indeed, from Eg. (11 001) we can see 
that at a fixed dark matter density profile the fifth-force 
enhancement of gravity by the factor (1-t-ei) yields an in¬ 
crease of the dark matter velocity dispersion and of the 
gas temperature by the same factor. The small devia¬ 
tions from this 2% value, which depend on mass, that 
appear in Fig. 1151 correspond to the small changes of the 
dark matter profile through the modification of the con¬ 
centration parameter shown in Fig. 1121 
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FIG. 16: Upper panel: temperature-luminosity relation for 
the K-mouflage models and the A-CDM reference, at z = 
0.048 (lower curves) and z = 0.451 (upper curves). The data 
points are taken from observations made by [39j (in green) and 
[33| (in brown), with clusters in the redshift range 0.048 < 
z < 0.451. Lower panel: relative deviation of the cluster 
temperature-luminosity relation from the A-CDM reference, 
at 2 = 0.048. 


At fixed temperature, the K-mouflage models give a 
slightly lower X-ray luminosity. This is because at fixed 
mass K-mouflage models give a higher temperature, as 
seen in Fig. 1151 Therefore they give a lower mass at 
fixed temperature. Since the X-ray luminosity scales as 
Lx ~ PsMVT, a lower mass implies a lower luminosity 
(disregarding the impact on p s ). As expected, we find 
percent deviations as for the M — T scaling law. 

Thus, as for the quantities studied in previous sec¬ 
tions, we obtain percent deviations from the A-CDM 
scaling laws. Unfortunately, this is probably too small 
to be used as a meaningful constraint on these modified- 
gravity scenarios, in view of the observational and theo¬ 
retical uncertainties. Therefore, it is unlikely that cluster 
scaling laws can provide competitive constraints on such 
modified-gravity models, that must also pass very tight 
Solar-System bounds and satisfy larger-scale cosmologi¬ 
cal constraints associated with the growth of large-scale 




FIG. 17: Upper panel: cluster temperature function for the 
K-mouflage models and the A-CDM reference, at 2 = 0.05. 
The data points are taken from observations made by [4(J 
from a sample of clusters with 2 ~ 0.05. Lower panel: rel¬ 
ative deviation of the cluster temperature function from the 
A-CDM reference at 2 = 0.05 

structures or the evolution of the Hubble expansion rate 
(e.g., constraints from BBN). 

F. Cluster Temperature function 

Neglecting the scatter of the mass-temperature rela¬ 
tion, by combining the halo mass function described in 
Sec. IIV El with the mass-temperature relation obtained in 
Sec. IV PI and Fig. [15] we obtain the X-ray cluster tem¬ 
perature function 

n{T) = n(M) ^F- (107) 

In Fig.[T7]we show the temperature functions computed 
for the K-mouflage models together with the A-CDM 
case, evaluated at 2 = 0.05 to compare them with the 
observations obtained by lioj ]. 

We obtain a reasonable agreement with observations. 
As is well known, this also means that the cluster tern- 
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perature is a rather robust quantity (as compared for in¬ 
stance with the X-ray luminosity) and that it is not nec¬ 
essary to build very sophisticated models to recover the 
right order of magnitude. As shown in the lower panel, 
we now obtain deviations for the cluster number counts 
that are of order unity: the K-mouflage models can pre¬ 
dict twice or three times more high-T clusters than the 
A-CDM reference. As we have seen, this difference is not 
due to deviations in the cluster scaling laws, that is, in 
the intracluster medium, which are quite small, but to 
the amplification of the high-mass tail of the halo mass 
function already shown in Fig. [TUI Therefore, this result 
should be rather robust as it is directly related to the 
faster growth of large-scale structures in the K-mouflage 
scenarios. 


G. Sunyaev-Zel’dovich effect 

An indirect method to infer the properties of the clus¬ 
ters is to use the Sunyaev-Zel’dovich effect (hereafter 
SZE) [13|. It occurs when photons from the CMB in¬ 
verse Thompson scatter in the intracluster medium. The 
measured CMB temperature is then distorted with an 
amplitude proportional to the so-called Compton param¬ 
eter (see e.g. [41| ), 

y = [ n e a T ——^dl, (108) 

J m e c z 

where n e is the electron number density, T g is the gas 
temperature, m e is the electron mass, ctt = 6.65 x 10~ 25 
cm 2 is the Thompson cross-section and d l denotes the in¬ 
tegration along the line of sight. Following a common ob¬ 
servational practice, by integrating over the angular area 
of the cluster, defined for instance by the radius i?sooc 
associated with the density contrast of 500 with respect 
to the critical density, we define the integral Compton 
parameter 

[ rRaaoc knT 

Tsooc = ydn = rj 2 (z) / ATrr 2 n e {r)a T - —fdr, 

J Jo m e c 

(109) 

where rd(z) is the angular distance of the cluster located 
at redshift z. 

In FigJTH] we show the relations M g 500c — Y 500c r^(z) 
and T500C—nocore - TsoocU^-) for the K-mouflage models 
and the A-CDM reference, and we compare them to the 
observations made by pl| . Again, we obtain a reason¬ 
able agreement with observations and a small deviation 
between the different scenarios. The agreement is bet¬ 
ter for the M g — Y relation than for the T — Y relation, 
but the latter shows a very large scatter and is probably 
contaminated by large observational errors. 

In any case, as in Sec. IV El it appears that deviations 
of cluster scaling laws associated with modified-gravity 
scenarios are too small as compared with observational 
error bars and theoretical uncertainties to be competi¬ 
tive. However, number counts, whether in terms of the 
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FIG. 18: Integrated Comptonization within Rsooc as a func¬ 
tion of the gas mass (upper panel) and gas temperature (lower 
panel) for the K-mouflage models and the A-CDM reference. 
These different models cannot be distinguished in these fig¬ 
ures. We show our results for 2 = 0.16 (lower curves in the 
upper panel and upper curves in the lower panel) and z = 1.45 
(upper curves in the upper panel and lower curves in the lower 
panel). The data points are measures from a sample of clus¬ 
ters in the range 0.16 < z < 1.45 [4l| . 


cluster temperature or SZ parameter Y, could provide 
useful constraints. 


H. Dynamical and weak lensing masses 

Finally, we briefly comment on the dynamical and 
weak lensing masses of clusters and massive halos. In the 
unscreened regime, which applies to clusters and larger 
scales, the dynamics of matter particles (dark matter 
and the gas) is governed by the metric potential $, as 
in Eos. (1151) and (15U) . This gravitational potential is re¬ 
lated to the matter density through the modified Poisson 
equation (l45l) and it is equal to the standard Newtonian 
potential (but with a time-dependent Newton constant) 
multiplied by the factor 1 + e \. From observations of the 
dynamics in clusters one would then measure the dynam- 
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FIG. 19: Relative deviation of the cluster correlation function 
from the A-CDM reference for the K-mouflage models. We 
consider halos of mass 10 15 h _1 AfQ (upper solid lines) and 
10 13 /i _1 M q (lower dashed lines). 


ical mass 

M dyn = [l + Cl (t)]^M, (110) 

yo 

assuming General Relativity (GR) gravity with today’s 
Newton’s constant. On the other hand, the weak lensing 
potential $ w i that governs the deflection of light rays by 
the perturbations of the metric is 

$ + vR 

$wi = —= 'M, (HI) 

where we used Eq. CD- Therefore, weak lensing observa¬ 
tions of clusters would give the weak lensing mass 


Mens = -^M, (112) 

y o 

and the ratio between the dynamical and weak lensing 
masses reads as 


Mlyn 

Mens 


1 + £l (t). 


(113) 


We have shown the factor ei in Fig. [5] Thus, we find that 
the dynamical mass is greater than the lensing mass by 
about 2 %. As explained above, this is set by the value 
of 2/3 2 , which is constrained to be of order 2% or below 
because of cosmological and Solar System constraints [HJ 
(BBN constraint and bound on the time dependence of 
G). Therefore, this result on the ratio Mlyn/Mens gives 
an upper bound for its deviation from GR for all realistic 
K-mouflage models. 


cluster correlation function. Following [42|, |43}, we write 
the halo bias as 

b{M) = l + IL-Zl. (H4) 

Ol 

More accurate fitting formulas have been proposed for 
A-CDM cosmologies [Hj], but they involve free param¬ 
eters that might vary for different modified-gravity sce¬ 
narios. Moreover, numerical simulations find that the 
spherical collapse model (11141) provides reasonably good 
predictions that can fare better than more sophisticated 
models for rare and massive halos (45^ . which we focus 
on here. Therefore, Ea. (11141) should be sufficient for our 
purposes and provide a simple estimate of the impact of 
K-nrouflage models. Note that because clusters are not 
screened, the reasoning that leads to Ea. (lll4D in the A- 
CDM cosmology remains valid for K-mouflage scenarios, 
as the only change is the time-dependent effective New¬ 
ton’s constant as it would be defined from Ea. dMl) . This 
enters the bias (11141) through the different values of v{M) 
and S L , that we compute from the spherical collapse de¬ 
scribed in Secs. IIII Dl and ITVEl 

Combining the halo bias (11141) with the matter correla¬ 
tion function £(;r) shown in Fig. llll we obtain the cluster 
correlation function £ c i(a") displayed in Fig. 1191 The com¬ 
parison with Fig. [XT] shows that the cluster correlation 
function is much less affected by K-mouflage than the 
matter correlation and it can actually be slightly lower 
than in the reference A-CDM cosmology. This is because 
the amplification of gravity, associated with the greater 
effective Newton’s constant (l+ei(t))G(t)/Go, merely ac¬ 
celerates the growth of large-scale structures. This am¬ 
plifies the matter density power spectrum and correlation 
function, as well as the large-mass tail of the halo mass 
function, as seen in Figs, fill and [TUI However, this same 
phenomenon also implies that, at a fixed mass M, mas¬ 
sive halos are less rare and have a smaller bias b(M) [in 
particular, i/(M) becomes smaller in Ea. (11141) ]. This ef¬ 
fect almost cancels the increase of the underlying matter 
density correlation function. Therefore, it appears that 
the matter correlation function, measured for instance 
from weak lensing observations or galaxy surveys (using 
typical halos with a bias of order unity that is not sig¬ 
nificantly changed by K-mouflage), is a better probe of 
such modified-gravity scenarios than the cluster correla¬ 
tion function (or more generally, the correlation of rare 
objects). 


VI. COMPARISON WITH OTHER 
MODIFIED-GRAVITY MODELS 


I. Cluster correlation function 

In the previous sections we considered the internal and 
integrated properties of clusters, as well as their abun¬ 
dance. Another probe of cosmology is provided by the 


A. Some other modified-gravity theories 

Before we conclude this study of K-mouflage models, it 
is interesting to compare our results with other modified- 
gravity models that have been investigated in the litera¬ 
ture. The main scenarios that have led to detailed analyt- 
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ical and numerical studies are the f(R) theories, scalar- 
field models such as Dilaton and Symmetron models, and 
Galileon models [ifil - lisl ]. 

The f(R) models [3, 49l-l53|. can be recast as scalar- 
field models with an Einstein-frame action of the form 
©, with a standard kinetic term, an exponential cou¬ 
pling function A(<£h and a scalar-field potential V(y>) 
[54j. The Dilaton [55H57j and Symmetron |H],[59| models 
are also scalar-tensor theories of this form, but with dif¬ 
ferent coupling functions A(tp) and potentials V(<p) (and 
standard kinetic terms). Finally, Galileon models [60j-[63| 
also involve a scalar field, with a non-standard kinetic 
term (the scalar-field Lagrangian contains higher-order 
terms in Vip and V 2 (/j), but there is no coupling func¬ 
tion A(ip) (i.e., = g^). (Of course, it is possible to 

build more complex models that combine these various 
ingredients.) 

These different scenarios show different nonlinear 
screening mechanisms that ensure convergence to GR 
in the Solar System, Chameleon 0 (for f(R) models), 
Damour-Polyakov S] (for Dilaton and Symmetron mod¬ 
els) and Vainshtein p (for Galileon models) mechanisms. 
The theories that are closest to K-mouflage scenarios are 
the Galileon models, as their screening mechanism also 
relies on the nonlinear derivative terms of the scalar- 
field Lagrangian; but they also involve V 2 </? instead of 
V<y9 only, which gives rise to different scaling exponents, 
for instance for the Vainshtein and K-nrouflage screening 
radii as a function of the mass M of compact objects. 


B. Einstein and Jordan frames 

A first important difference between the K-mouflage 
scenario and these other modified-gravity models is the 
distinction between the Einstein and Jordan frames. As 
recalled above, this distinction does not apply to the 
(simplest) Galileon models, but the f(R), Dilaton and 
Symmetron models also naturally give rise to distinct 
Einstein and Jordan frames. However, it turns out that 
in these scenarios the coupling function is constrained to 
remain very close to unity. Thus, \A — 1| < 10~ 4 for 
f(R) models, because the “mass” of the scalar field must 
be sufficiently large, m > 10 3 H/c, to ensure an efficient 
screening of the fifth force by a Chameleon mechanism in 
the Solar System. For Dilaton and Symmetron models 
we have \A — 1| < 10 -6 as the coupling strength /3 must 
vanish sufficiently fast in high-density regions to screen 
the fifth force through a Damour-Polyakov mechanism. 
This means that, in terms of background quantities (e.g., 
the Hubble expansion rate and the scale factor) one can 
identify the Einstein and Jordan frames, which also be¬ 
come almost identical to the A-CDM reference. However, 
at the level of the metric perturbations <f> and 4', this is no 
longer the case and the Einstein and Jordan gravitational 
potentials differ by terms set by <L4, and the dynamics of 
perturbations deviate from the A-CDM reference because 
of the fifth force. 


In the K-mouflage case, this identification already 
breaks down at the background level. Indeed, \A— 1| can 
reach values of the order of a few percents (see Fig. © 
while being consistent with Solar System and cosmologi¬ 
cal constraints [lj|] ■ For the same reason, the background 
dynamics (in both Einstein and Jordan frames) show 
percent deviations from the A-CDM reference. There¬ 
fore, we must pay attention to the distinction between 
Einstein and Jordan frames already at the background 
level. In particular, in this paper, as we study clusters of 
galaxies that involve atomic or radiative processes (both 
for the definition of their redshift, from atomic lines, and 
for their properties such as X-ray emission), the Jordan 
frame is the one that is more directly connected to obser¬ 
vations and we work in this frame. Another advantage 
of the Jordan frame is that the equations of motion take 
their usual form, in particular matter is conserved, which 
permits a clear and simple physical interpretation, and 
only gravity is modified. In contrast, in the Einstein 
frame gravity takes a standard form but the equations 
of motion are modified and the matter density is usually 
not conserved. 


C. Scale dependence and screening regime 

A second important difference between the K-mouflage 
scenario and some other modified-gravity models is that 
the deviations from the A-CDM reference are scale inde¬ 
pendent on perturbative scales (from cluster scales to the 
horizon). This is most easily seen from the fact that the 
factor ei (t) that enters the evolution equation (1551) of the 
linear growing mode only depends on time, so that the 
linear growing mode D + (t) remains scale independent as 
in the A-CDM cosmology. This is due to the fact that 
in the scalar-field Lagrangian ([6]) we focused on the non¬ 
standard kinetic term and neglected a possible potential 
term V(<p). Of course, in the highly nonlinear regime a 
new scale dependence appears, as the fluctuations of the 
scalar field themselves become nonlinear and give rise 
to the K-mouflage screening mechanism, which ensures 
the convergence to GR in the Solar System. However, it 
happens that the screening transition appears at galaxy 
scales, so that clusters remain unscreened and fully feel 
the effect of the fifth force. 

The same behavior is obtained in the Galileon mod¬ 
els, where linear scales below the horizon show scale- 
independent growing modes and the Vainshtein screening 
mechanism applies to cluster scales and below |64j, HU ; 
but in the K-mouflage case the nonlinear screening only 
applies to galaxy scales and below, as clusters remain un¬ 
screened. In contrast, in f(R) and Dilaton/Symmetron 
models, there is a characteristic scale dependence, as we 
recover GR both on very large scales x lh -1 Mpc (be¬ 
cause of the finite mass of the scalar field) and on very 
small scales x <C lh -1 Mpc (because of nonlinear screen¬ 
ing mechanisms, here Chameleon or Damour-Polyakov 
mechanisms S 0). Then, the linear growing mode 
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D + (k,t) shows a clear scale dependence on quasilinear 
scales and nonlinear screening effects also add a further 
scale dependence around x ~ l/i -1 Mpc [gfif]. Thus, in 
these models clusters probe a scale-dependent regime 
and the transition between the unscreened and screened 
regimes. 

Therefore, clusters of galaxies are especially well-suited 
probes of K-mouflage scenarios because they are un¬ 
screened (hence they feel the full amplitude of the fifth 
force). Moreover, the modification of gravity is still scale 
independent on these scales so that cluster properties 
should not be too difficult to model (the same modeliza- 
tions should apply equally well to the A-CDM and K- 
mouflage cosmologies). 


D. Clusters 

1. Cluster profles 

The effect of the fifth force on the cluster matter and 
gas profiles within the context of modified-gravity scenar¬ 
ios with chameleon mechanisms [mostly for f(R) models] 
has been investigated in [67l - l70| . 

As recalled above, the fifth force effect is somewhat 
different between f(R), Dilaton and Symmetron models, 
and K-mouflage scenarios. In the former cases clusters 
are typically at the transition between the screened and 
unscreened regimes. Then, massive clusters are screened 
by the Chameleon or Damour-Polyakov mechanisms (the 
deviations from GR being most efficiently suppressed in 
the Symmetron models) while for low-mass clusters only 
a small core is efficiently screened. In particular, for f(R) 
and Symmetron models the amplification of gravity is lo¬ 
calized in the outskirts of massive clusters. This gives a 
distinct scale dependence of the modified-gravity effect 
in these theories 0lll, but the efficient screening also 
decreases the overall deviation from GR. See for instance 
jUJ for a detailed analysis and comparison of these dif¬ 
ferent models and (zj, l73l - l76j for numerical simulations 
of various models. In Galileon scenarios, depending on 
the model, clusters may be fully screened (and their pro¬ 
files are similar to those of a quintessence model with no 
fifth force and the same expansion history) or only partly 
screened (which gives rise to complex effects) |65[. In 
the K-mouflage case, clusters still are in the unscreened 
regime. Therefore, there is no characteristic scale depen¬ 
dence that can be used to distinguish them from A-CDM 
cosmologies, but the amplitude of the smooth deviation 
is greater (as compared with screened scenarios). 

There have been no specific simulations of K-mouflage 
models so far, but we can recall here some results from 
simulations of other modified-gravity models. 

Ref. [z3 develops hydrodynamic N-body simulations 
to investigate the impact on dark matter and gas profiles 
of f(R) and Symmetron scenarios. In agreement with the 
discussion above, the authors find that the dark matter 
density is increased as compared with the A-CDM refer¬ 


ence in the outskirts of massive halos. This is because the 
fifth force applies to outer radii, which are unscreened, 
and this also yields a greater velocity to the particles, 
which cannot cluster as strongly within inner radii. They 
also find a lower deviation from A-CDM for the gas den¬ 
sity profile than for the dark matter density profile. They 
note that this may be due to delays in the collapse of the 
dark matter and the gas, with the screening of the halos 
taking place after dark matter collapse and before gas col¬ 
lapse. As explained above, for K-mouflage scenarios we 
do not expect such a localized enhancement in the matter 
densities and different behaviors for the dark matter and 
the gas, because clusters are unscreened and the dynam¬ 
ics remains similar to the A-CDM cosmology, as illus¬ 
trated by the explicit symmetry described in Sec. ITTTcl 
In fact, in Fig. [IT] we find that within our very simple 
model the deviation from A-CDM is slightly greater for 
the gas than for the dark matter [because of the small 
change in the concentration parameter and the higher 
sensitivity for the gas that arises from the exponential 
equilibrium distribution (1981) ]. 

We can note that N-body simulations of f(R) models 
also find that in the case f Ro = —10~ 4 , where Chameleon 
screening is not efficient, deviations from A-CDM are 
smooth and the velocity dispersion and gas tempera¬ 
ture T g are about 4/3 times the GR value [H, [zl[, due 
to the 4/3 increase of the effective Newton’s constant. In 
our case, this corresponds to the (1 + ei) enhancement 
in Eci. dlOOD . but with e\ ~ 2% instead of 1/3. These 
simulations also find that dark matter halos remain well 
described by NFW profiles [77]. These results suggest 
that our approach for the density profiles and the gas 
temperature described in Secs. IV Al and IV Dl should fare 
reasonably well. 

For Chameleon scenarios, Ref. [67] considers the effect 
of the fifth force on the cluster gas profile, which be¬ 
comes more compact for a given dark matter profile [in 
our case, we also include the effect on the dark matter 
profile, through the modification of the concentration pa¬ 
rameter shown in Fig. 1121 and we assume an isothermal 
gas (i.e., 7 = 1) instead of a polytropic equation of state 
p g oc with 7 ~ 1.2]. However, they find that observa¬ 
tional error bars are too large to give useful constraints 
on f(R) models. We reach the same conclusion for the 
K-mouflage models studied in this paper. 

Refs. [H,|7g combine X-ray observations (which probe 
the temperature and electron number density profiles) 
and weak lensing signals (which probe the total mat¬ 
ter profile) to constrain deviations from General Relativ¬ 
ity. Indeed, while the gas profile is sensitive to the fifth 
force the lensing deflection of light rays remains the same 
as in GR. This allows them to derive the upper bound 
\f Ro \ < 6 x 10 -5 . The same behavior applies to the K- 
mouflage scenario. In this paper we found a few percent 
deviations from GR, which should apply to all realistic 
K-mouflage models that satisfy cosmological and Solar 
System constraints [as this is due to the constraint on 
the coupling /?, independently of the details of the kinetic 
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function K (%)]. We leave a more general MCMC analy¬ 
sis of K-mouflage scenarios, combining different probes, 
to future works. 


2. Cluster scaling relations 

The impact of f(R ) gravity on the cluster scaling re¬ 
lations has been studied in Ref. [&§] with numerical sim¬ 
ulations. Again, we find similar behaviors for the K- 
mouflage models as for the f(R) model with fn 0 = 
—10“ 4 , where clusters are not screened. In particular, 
the dark matter velocity dispersion and gas temperature 
are increased, at fixed mass, and the X-ray luminosity is 
decreased, at fixed temperature, as compared with the A- 
CDM reference. However, whereas in f(R) theories the 
deviations from A-CDM can reach a factor 1 /3 in the un¬ 
screened regime, and deviations of order unity can also 
be expected in Dilaton or Symmetron models, realistic 
K-mouflage models can only deviate by a few percent at 
most. Indeed, the magnitude of these deviations is set 
by the factor ei, which itself is set by 2/3 2 (at z = 0) 
and the coupling strength must satisfy (3 < 0.1 because 
of cosmological and Solar System constraints [l3 |. 

3. Cluster lensing; Dynamical and leasing masses 

In f(R), Dilaton and Symmetron scenarios the weak 
lensing potential $ w i that governs the deflection of light 
rays, given by $ w i = (<h + /2, is equal to the Newto¬ 

nian potential Tn [using \Jr 0 \ <C 1 in the case of f(R) 
theories] [78]. Then, for weak lensing observations the 
only difference from the A-CDM scenario comes from the 
different evolution of the matter density fields. In con¬ 
trast, the motions of matter particles or of small halos 
(e.g., satellite halos or small galaxies), which fall towards 
massive clusters, feel the fifth force. This gives rise to 
different lensing and dynamical masses. This has been 
investigated in semianalytical models and N-body simu¬ 
lations [xl, [z3 and used in [H, Q in combination with 
X-ray observations, for f{R) theories, and in jgflj, using 
the galaxy infall kinematics onto massive clusters, for 
f(R) and Galileon models. 

In the K-mouflage scenario, we again have <I> w i = 'I'n, 
as the factor 6 A/A cancels from the sum over the two 
metric potentials, see Eq. m • However, in contrast with 
the former theories the effective Newton constant [which 
enters the Newtonian potential flU] now depends on 
time, as in Eg. (1771) . (This effect does not appear in the 
other scenarios because they have A very close to unity, 
within 1CT 4 or better.) On the other hand, this time- 
dependent prefactor cancels from the ratio between the 
dynamical and lensing potentials or masses, see Eci. (11131) . 
However, whereas in f(R) theories this ratio can again 
deviate from unity by a factor 1/3 in the unscreened 
regime, and deviations of order unity can also be reached 
in Dilaton or Symmetron models, in realistic K-mouflage 


models this ratio can only deviate from unity by 2% at 
most because of observational constraint on the scalar- 
field coupling, /3 < 0.1. Therefore, K-mouflage models 
cannot significantly decrease the tension between mea¬ 
sures of X-ray and lensing clusters masses. 

4. Cluster number counts 

In most modified-gravity scenarios the growth of large- 
scale perturbations differs from the GR evolution. This 
typically leads to a new scale dependence of the linear 
growth rates (e.g., in /(!?), Dilaton and Symmetron mod¬ 
els), as one goes from the very large scales (beyond the 
Compton wavelength of the scalar field) where GR is re¬ 
covered to the quasi-linear scales where the fifth force 
is unscreened and gives its maximum amplification of 
the gravitational interaction (at smaller scales nonlin¬ 
ear screening leads again to a convergence to GR). This 
amplification typically yields a faster growth of matter 
density perturbations on scales lft -1 Mpc to 10/i~ 1 Mpc, 
whence a greater abundance of massive halos and clus¬ 
ters as compared with the A-CDM cosmology, see Refs. 
[65 [73l IsTUsI ] for numerical studies of various models. 
As explained above, a similar enhancement is found in 
the K-mouflage scenarios, with the important difference 
that all clusters are unscreened and that the modification 
of gravity extends up to the horizon (so that the linear 
modes grow faster than in A-CDM but remain scale in¬ 
dependent). In Galileon scenarios, the screening mecha¬ 
nism has a strong impact and, depending on the models, 
the tail of the halo mass function can be either increased 
or decreased, as compared with a quintessence scenario 
with the same expansion history [65[. 

The abundance of massive clusters has been used 
within f(R) theories to constrain /r 0 [85l - [87l . In com¬ 
bination with CMB, BAO and SNela observations, one 
can obtain an upper bound \fn 0 \ < 1.7 x 10 -5 86. 87], 
but most of the constraint comes from the cluster data. 
We can expect that for K-nrouflage models similar re¬ 
sults should be obtained, especially since clusters are un¬ 
screened so that their abundance should provide useful 
constraints. In this paper we presented the physics of K- 
mouflage scenarios on cluster scales, highlighting the dif¬ 
ference between the Einstein and Jordan frames (which 
can be neglected in most other scenarios) and investigat¬ 
ing both the modified growth of structures and halo mass 
functions and the modified cluster scaling laws. We leave 
a detailed MCMC analysis of the K-mouflage parameter 
space to future works. 


5. Cluster correlation function 

In modified-gravity scenarios, the correlation function 
and the power spectrum of the matter density field are 
often increased as the growth of large-scale structures 
is amplified by the fifth force. This also enhances the 
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large-mass tail of the halo mass function and decreases 
the bias of massive halos as they become less rare. In the 
K-mouflage models that we considered in this paper, this 
decrease of the cluster bias mostly cancels the increase of 
the underlying matter density correlation, and the clus¬ 
ter correlation function is much closer to the A-CDM 
reference than the matter correlation itself. Therefore, 
the correlation of the matter density field, which can be 
measured from weak lensing observations for instance, is 
a better probe of modified cosmology than the clustering 
of massive halos. The same effects, with a similar com¬ 
pensation between the smaller bias and the higher mat¬ 
ter correlation (but the compensation may only be par¬ 
tial, depending on the model) appear in other modified- 
gravity scenarios, see for instance Ref.|§8} for f(R) the¬ 
ories. 


E. Other tests 

On slightly larger scales than those probed by clus¬ 
ters of galaxies, modified-gravity models have also been 
tested from galaxy surveys, using redshift-space distor¬ 
tions of the galaxy power spectrum [89} , the clustering of 
Luminous Red Galaxies 63, or the shape of the galaxy 
correlation function itself [90}. Similar studies could be 
performed for K-mouflage models, as they can lead to 
20% deviations for the matter power spectrum and cor¬ 
relation function, as seen in Fig. [TT| This will be inves¬ 
tigated in future works. 

In between the cluster and cosmological scales (e.g., 
the formation of large-scale structures and the back¬ 
ground dynamics) and the Solar System scales (where we 
must recover GR up to a very high accuracy), modified- 
gravity theories have also been tested on intermediate 
galaxy scales. In particular, within Chameleon scenar¬ 
ios, low-mass galaxies should be unscreened so that the 
rotation curve of their diffuse gas component probes the 
fifth force, while their stars, being compact high-density 
objects, are screened and move as in GR [H}. This can 
provide constraints as tight as | /_r 0 | < 10 -6 for f(R) 
theories [92}. The same behaviors apply to K-mouflage 
models and we can expect that such tests could also pro¬ 
vide useful constraints on these scenarios. We leave this 
task to future studies. 

Numerical simulations have also been used to inves¬ 
tigate the impact of f(R) theories on Lyman-a forest 
observations (the transmitted flux fraction and the flux 
power spectrum) [93]. They find changes that are too 
small as compared with observational error bars and 
do not provide competitive constraints. Although this 
study does not distinguish the impact of screening ef¬ 
fects, screening can be expected to be rather inefficient 
for such moderate density fluctuations. Therefore, it is 
likely that similar conclusions would be reached for K- 
mouflage models, but we leave a detailed study for future 
works. 


VII. CONCLUSIONS 

We have extended previous works on K-mouflage mod¬ 
els by deriving the equations of motion in both the Ein¬ 
stein and Jordan frames for a fluid with pressure, and 
next focusing on the usual case where the pressure arises 
from small-scale nonlinear processes. In contrast with 
many modified-gravity scenarios, the Einstein and Jor¬ 
dan frames already differ by a few percent at the back¬ 
ground level, for K-mouflage models that are consistent 
with cosmological and Solar System constraints. There¬ 
fore, one must take into account these deviations and use 
the correct quantities when comparing with observations. 

We focused on the Jordan frame, which is better suited 
to cosmological probes that involve atomic processes, 
such as X-ray clusters. We show that even though K- 
mouflage models only differ from the A-CDM reference 
by a few percent at the background level, the linear grow¬ 
ing modes can deviate by 10% and the matter density 
power spectrum and correlation function by 30%. The 
tail of the halo mass function is enhanced by factors of 
order unity for M > 1O 14 /i _ 1 M 0 at z = 0. 

Galaxy clusters are not screened by the K-mouflage 
mechanism, contrary to what happens for chameleon 
models like /(i?) in the large curvature limit or Galilcon 
models subject to the Vainshtein screening. For this rea¬ 
son, we investigate the effects of the K-mouflage mod¬ 
ification of gravity on clusters. We take into account 
the impact of the fifth force mediated by the K-mouflage 
scalar field on both the dark matter and gas profiles, 
through the modifications to the growth of large-scale 
structures and to the hydrostatic equilibrium. We find 
that K-mouflage makes clusters slightly more compact. 
Similarly, the gas temperature and the X-ray luminos¬ 
ity differ from their A-CDM counterparts by a few per¬ 
cent, an order of magnitude which follows directly from 
the constraints on K-mouflage (especially on the coupling 
strength /3) due to Solar System tests. The only deviation 
of noticeable order appears in the cluster temperature 
function, as the number of clusters is larger than in the 
A-CDM scenario for K-mouflage models, because of the 
faster growth of large-scale structures. This appears to 
be large enough that one can hope that this will be within 
the reach of the future large-scale surveys. On the other 
hand, the cluster correlation function only deviates by a 
few percent from the A-CDM case because the increase 
of the underlying matter density correlation function is 
compensated by the lower cluster bias, as massive ha¬ 
los become less rare at fixed mass, due to the enhanced 
structure formation. 

In this paper we only considered two kinetic functions 
K(x) to illustrate the K-mouflage scenarios and to esti¬ 
mate the amplitude of the effects that can be reached by 
realistic models, which are consistent with cosmological 
and Solar System constraints. We leave to future works a 
more detailed MCMC analysis of K-mouflage scenarios, 
which would provide the parameter space of K-mouflage 
models that is allowed by cluster observations and com- 
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binations with other observational probes. 
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Appendix A: Equations of motion in the Einstein 
frame 

In this appendix we derive the equations of motion of 
the scalar field and of the matter component in the Ein¬ 
stein frame, for a cosmological fluid with a nonzero pres¬ 
sure. The derivation is similar to the one presented in 
previous papers }7ii, H|, where we studied the background 
cosmological dynamics and the formation of large-scale 
cosmological structures, but with the addition of the 
pressure terms. 

1. Energy-momentum tensors and equations of 
motion 


Here we used the fact that the trace of the matter energy- 
momentum tensor is TT = — p + 3 p, from Eq. m, while 
it is zero for the radiation component as p( r ) = P( r )/3. 
Combining with Eq. (IA2l) . we find for the scalar field the 
“non-conservation” equation 

V^ ) = (p-3p)V,lnA (A4) 

The matter energy-momentum tensor is conserved in 
the Jordan frame, 

= 0, (A5) 

which gives in the Einstein frame the non-conservation 
equation 

= — (p — 3p) Vj, In A. (A6) 

On the other hand, the radiation energy-momentum ten¬ 
sor is conserved in both frames, 

= o, V„f; (r) = 0. (A7) 

The sum of all energy-momentum tensors is also con¬ 
served in both frames, and + Tj^ + T^, = 

V^ + ^ (r) +T^ } ]=0. 

Finally, in the Einstein frame the Einstein equations 
take the usual form, Mp\G£ =T^+ Tj^ + T^. 


We consider three components of the energy den¬ 
sity of the Universe, a matter fluid with nonzero pres¬ 
sure, radiation, and the scalar field. The Einstein-frame 
and Jordan-frame matter energy-momentum tensors are 
given by Ea. ITOl) and they satisfy the relations (fill) . As¬ 
suming a perfect fluid, the matter energy-momentum ten¬ 
sor writes as 

f£ = (p + p)u li u v +p5jZ, (Al) 

where u^ is the velocity 4-vector, normalized such that 
u^Ufi = — 1, and p and p are the Einstein-frame density 
and pressure, which are related to their Jordan-frame 
counterparts by Eq. m- 

For the radiation component, we neglect perturbations 
and only consider the contribution to the background, 
with the mean density and pressure p( r ) and p( r ) = p( r )/3. 
Their Jordan-frame counterparts are again given as in 
Eq.(®. 

The Einstein-frame energy-momentum tensor of the 
scalar field reads as 


2. Background dynamics 


At the level of the cosmological background, the Fried¬ 
mann equations take the usual form in the Einstein 
frame, 

ZMp\H 2 = p + p( r ) + Ptp, (A8) 

dH 

— 2Mpj— = p + p + p( r ) + P( r ) + pip + P v , (A9) 

where p v and p v are the background scalar-field energy 
density and pressure (in the Einstein frame), given by 
Eq.tflfl). 

The Klein-Gordon equation (IA3D gives 


d_ 
d t 


s 3 Sa" 

d£ 


-a J (p- 3p) 


d In A 


dtp 


and the scalar-field energy density satisfies 


(A10) 


U) = J=§^ = *"V M ¥>V„¥> + M A Kg^. (A2) 

The Klein-Gordon equation that governs the dynamics 
of the scalar held <p is obtained from the variation of the 
action m with respect to <p. This yields 


=-3H(p v +p ip )-(p-3p)^^. (All) 
d t at 

The non-conservation equation (IA6I) gives for the mat¬ 
ter density the evolution equation 


Wtp K' = (p - 3 p) 


d In A 

dtp 


(A3) 


dp 
d t 


~ - —. „ - d In A 

-3 H(p + p) + (p- 3p)- 


d t 


(A12) 
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In particular, we have d(p + p^) / dt = — 3H(p + p + p v + 
p v ). When the pressure is zero, we can define a conserved 
density p by 


p = p/A, p=p/A. (A13) 

Indeed, substituting into Ea. (IA12l) gives 

W = (A14) 

Thus, if p = 0 we obtain a conservation of the stan¬ 
dard form in the Einstein frame, dp/dt = —3Hp. How¬ 
ever, if p ^ 0 it is no longer possible to cancel the non¬ 
conservation term of Ea. (IA12l) by such a simple redefini¬ 
tion of the density. 

The background radiation density obeys the usual con¬ 
servation equation, 

= -3#(/5(r) + m ) = — 4^( r) , (A15) 

in agreement with the second Eq. (1A71) . 


3. Perturbations 


Contracting with u v and using vyu u = —1, we get 

p + 3h(p+p) = (p- 3p)^. (A21) 

It is easy to see that p = a~ l [d T p + (v • V)p], where 
V = d/dx is the standard 3D spatial gradient, and 3 h = 
d _1 [3"H + (V ■ v)], where H = dlnd/dr is the conformal 
expansion rate in the Einstein frame. Therefore, this is 
explicitly 


+ (v • V)/5+ (3U + V • v)(p +p) 


= (P~ 3p) 


a in a 

~d7~ 


+ (v- V)lnA 


(A22) 


Next, the non-conservation equation (IA19I) can be sim¬ 
plified by subtracting Ea. (IA21D multiplied by u v . This 
leads to 


V„p + u u p p — 3p V„A + u u A 


P + P 


P+P 


A 


(A23) 


This is the generalized geodesic equation. Specializing to 
p, = i, we get the Euler equation of K-mouflage hydrody¬ 
namics, 


In the conformal Newtonian gauge the Einstein-frame 
metric can be written as 

ds 2 = d 2 [—(1 + 2<E>)dr 2 + (1 - 24')dx 2 ], (A16) 

where r = J dt/a is the conformal time, and $ and 
are the two metric gravitational potentials. Throughout 
this paper, we consider the non-relativistic and weak- 
gravitational-fields regime, with <C 1, T <C 1, v 2 <C 1 
(where v = dx/dr is the peculiar velocity of the par¬ 
ticles), and we expand up to first order in {<E>, T, u 2 }. 
Then, in the final equations, we only keep zeroth-order 
terms, 1 + (h ~ 1 and 1 + v 2 ~ 1, except when the po¬ 
tentials or the velocity arise with a gradient operator, as 
in Eas. (|A24|) and (1A27I) . In particular, we have for the 
matter velocity 4-vector, 

u M = a _1 (l — $ + u 2 /2, id), = —d(l-t-$ + w 2 /2,-Uj), 

(A17) 

where we denote 

Vi = v l = — (A18) 


the peculiar velocity. 

The non-conservation equation (IA6I) gives 

(p + p)u v + 3h(p + p)u v + (p + p)u v + = 

— (p — 3p) V u (In A), (A19) 

where we have introduced 

p = u^V^p, u v = pUv, 3h = (A20) 


civ 
d t 


+ (v- V)v + 



1 dp 
p + p dr 


p — 3p <9 In A 
p + p dr 


v = 


-Vi - (A24) 

p+p p+p 


From Eci. (lA24|) we can see that (v • V) In A ~ d T v 2 + 
u 2 (V • v), hence this term can be neglected in the conti¬ 
nuity equation (IA22I) in the non-relativistic limit v 2 -C 1. 
This simplifies as 


^ + (v.V)p+(3-H+V-v)(p+p) = (p- 3 p)^—. (A25) 
The Klein-Gordon equation (1A3I) writes as 


l_d_ 
a 3 dt 



4 V • (Vp K') 
a z 


~(P~ 3p) 


din A 

dp 

(A26) 


The (0,0) component of Einstein equations, M^Gq = 
T 0 ° + T 0 %), gives 


’h = ’I'm with -^V 2 'k N = —l— (Sp + Sp^), (A27) 

a 2 2Mp, 

where Sp = p — p and Sp^ = p v — p v , with p v = 
—M^K + K r (dp/dt) 2 . Here we denoted Tn the usual 
Newtonian potential. The ( i,j ) components of the Ein¬ 
stein equations give (focusing on the part that is not 
proportional to ). 

M 2 ~ - K' 

-r/p-didj (\fr - $) = (p + p)viVj + -r^dipdjp. (A28) 
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From Eci. (IA27|) we can see that on the left-hand side in 
Ea. (lA28[) we have Mp l a~ 2 didj'i l ~ Sp ~ p. Therefore, 
the first term on the right-hand side, (p + p)uiVj, is neg¬ 
ligible as v 2 <C 1. Next, from the Klein-Gordon equation 
(IA26I) we obtain, on scales that are much smaller than 
the horizon, 

where k is the typical comoving wave number of interest. 
Then, the second term in the right-hand side in Ea. (IA28|) 
is of order 

K' . Sp /3 2 a 2 H 2 , Aon . 

-dt<pdj<p ~ Spj— 7 -p- « Sp, (A30) 

which is again negligible compared to M 2 - >x a^ 2 () l d : j x l>. 
Therefore, the Einstein equations (IA28I) give 

= T = 'I'n, (A31) 

within the approximations that we use in this paper. 

To close the system formed by the equations of mo¬ 
tion obtained above, we must specify the pressure, for 
instance through an (effective) equation of state such as 
p = wp with some parameter w. 


4 . Pressure due to small-scale nonlinear physics 

In the previous sections, we derived the equations of 
motion for a cosmological fluid with a nonzero pressure, 
in the non-relativistic limit v 2 <C 1 for the mean fluid 
velocity and in the weak field regime $jj€l. We made 
no approximation for the pressure and the equations of 
motion also apply to fluids with a pressure of the same 
order as the density, such as p = wp where w is a param¬ 
eter of order unity. However, in the usual CDM context, 
the pressure is negligible on cosmological scales and it is 
built on small scales by nonlinear processes, such as the 
collapse of gas clouds that generate shocks or the virial- 
ization of dark matter halos (which generate an effective 
pressure through the velocity dispersion of the particles). 
Then, the pressure is of the order p ~ pc 2 , where c s is the 
speed of sound or the velocity dispersion, and c 2 ~ Tn 
because it is generated by the gravitational collapse (for 
instance, if we have hydrostatic equilibrium we typically 
have V# N ~ Vp/p as the pressure balances gravity). 

Then, in the regime p/p ~ ’Fn ~ v 2 <C 1, the back¬ 
ground pressure is zero, p = 0, and we recover the cos¬ 
mological dynamics studied in Q for a pressureless fluid. 
Thus, the Friedmann equations read as Eqs. m-m and 
the matter and radiation densities evolve as in Eci. dlUl) . 
The Klein-Gordon equation becomes as in Eq. m- 

For the perturbations, the continuity and Euler equa¬ 
tions (IA25I) and (IA24I) simplify as 

+ V ■ (pv) + 3 Up = (A32) 


and 

dv (~ <91nA\ ~ Vp 

—+ (v-V)v+ in + ^ j v = -V^N + lnA)-—, 

(A33) 

while the Poisson equation remains identical to Eq. (1A27I) 
and the Klein-Gordon equation (IA26D becomes 


J_d_ 

a 3 dt 



4v-(VpK') = -p^. (A34) 
a z dtp 


Therefore, in this regime the only effect of the pressure 
is to add the usual pressure term in the Euler equa¬ 
tion, without mixed terms involving the coupling func¬ 
tion A(ip). 


5. Sub-horizon regime 


To simplify the Einstein equations (IA28I) we already 
used the small-scale limit, k/aH 1, which corresponds 
to scales that are much below the Hubble scale fu = 1 /H. 
This is the regime that is relevant for the formation 
of cosmological large-scale structures, such as clusters 
of galaxies. Then, the continuity and Euler equations 
(IA32I) - (|A33I) and the Poisson equation (IA27I) can be fur¬ 
ther simplified. Indeed, as in Ea. (IA29l) . we obtain the 
estimates 


k/aH 1 
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Mpi 

fSSp p 2 a 2 H 2 Sp 


P a 2 H 2 Sp ^ Sp 

~K' k 2 J ^ TP 
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K' k 2 


Sp 

< -S-, 
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S 'X ^ - 
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(VM 2 


P 2 ~a 2 H 2 (Sp) 2 (Sp) 2 

~o 1 


2M a cl 2 

K' 2 k 2 

p 2 

P 2 ct 2 H 2 

(i + *£) 

< 1. 

K' k 2 

\ P ) 


(A35) 

(A36) 

(A37) 

(A38) 


Then, in the continuity and Euler equations (IA32I) - (IA33D . 
we can write dhiA/dr ~ din A/dr, which leads to 
Eqs. and (12T1) . In the Poisson equation (IA27I) , we can 
neglect Sp v , which gives Eq. lB^l) . In the Klein-Gordon 
equation (IA34I) , we can neglect the fluctuations of A and 
only keep the spatial gradients. This leads to Ea. (E^l) . 
which also corresponds to the quasi-static approximation. 


6. Formation of large-scale structures 

Introducing the Einstein-frame matter density con¬ 
trast, 

6 = Sp/p, (A39) 


the continuity equation (12T71) also writes as 
dS 


dr 


+ V • [(1 + (S)v] = 0. 


(A40) 
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Thus, in terms of the density contrast we recover the 
usual continuity equation, without any A-term left. This 
is related to the fact that 5 = 6 , where 5 = Sp/p is the 
conserved matter density introduced in Ea. (IA13|) . within 
our set of approximations (p <C p and A ~ A, so that the 
factor A cancels out in the ratio 5p/p). 

On perturbative scales, we set the pressure term to 
zero, as in standard perturbation theory, because it is 
generated by non-perturbative effects such as shell cross¬ 
ing and virialization (shocks). Then, the formation of 
large-scale structures can be tackled through a perturba¬ 
tive approach, as in the usual A-CDM case. Introducing 
the two-component vector ip, 

($0 - (-(V-v)/(do/dt))’ (A41) 

equations (1A40I) and (12T1) read in Fourier space as 

- V>2 = /dkidk 2 <5D(ki+k 2 -k)d(ki,k 2 ) 

ama J 

x V> 2 (ki)V>i(k 2 ), (A42) 

Mk - ^ m(1 + + ( 2 ++e ' 2 ) ^ = 

dk 2 6 d (k, + k 2 - k) $ (k,, k 2 )ip 2 (kr )^2 (k 2 ),(A43) 


with 


d(ki,k 2 ) 


(ki + k 2 ) -kj 

kf 


/3(ki,k 2 ) 


|ki + k 2 | 2 (kj-k 2 ) 

2k 2 k 2 


(A44) 

The two differences from the equations of motion ob¬ 
tained in the A-CDM cosmology are the two time- 
dependent factors defined by 



h{i) = 


din A 
d In a 


P dy 

M P1 d In a 


(A45) 


In Eci. (lA43l) the factor f2 m (l + ei) can also be written as 
flm(l-t-ei), where Cl m is the cosmological parameter asso¬ 
ciated with the conserved density p defined in Ea. (IA13l) . 

[il'm = P /Pcrit = f^m/A, with Pcrit — 3.1/j , | / / the 

Einstein-frame critical density], and = A(1 + ey) — 1 = 
A — 1 + 2 A/3 2 /AT 

On large scales or at early times, we can linearize the 
equations of motion (IA42I) - (IA43I) . This gives for the lin¬ 
ear growing and decaying modes D±(t) the evolution 
equation 


d 2 D 

d(lna) 2 


1 d H 
~FP~dt 


£2 


J T - —^^m(l + ei)£ ) — 0. 

d in a 2 


(A46) 
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FIG. 20: Relative deviation, H/H — 1, of the Einstein-frame 
Hubble rate from the Jordan-frame Hubble rate, as a function 
of the Jordan-frame redshift. 



z 


FIG. 21: Deviation of the Einstein-frame cosmological pa¬ 
rameters from their Jordan-frame counterparts, Cl — D, as a 
function of the Jordan-frame redshift. We show the matter 
density parameters (solid lines) and the dark-energy density 
parameters (dashed lines). In the Einstein frame we consider 
the effective matter and dark-energy densities as given by 
Eas. (IB2l) and (IB3I) . 


Appendix B: Comparison of Einstein-frame and 
Jordan-frame backgrounds 

In this appendix we compare the Einstein-frame back¬ 
ground quantities with their Jordan-frame counterparts. 
We show in Fig. [221 the relative deviation of the Einstein- 
frame Hubble rate from the Jordan-frame Hubble rate. 
From Ea. (l24l) this is given by 


^ = A(1 - e 2 ) - 1. (Bl) 

We can see that at low redshift the Einstein-frame Hubble 
rate is typically lower than its Jordan-frame counterpart, 
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at a fixed Jordan-frame redshift. The comparison with 
Fig. [2] shows that at z < 4 the deviation between these 
two expansion rates is greater than the deviation between 
the Jordan-frame expansion and its A-CDM reference. 
As expected, this deviation is again of the order of a few 
percent, set by the value of /3 2 . This clearly shows that 
for the K-mouflage scenario one cannot treat both frames 
as approximately identical, contrary to what happens in 
many modified-gravity theories such as f(R) models or 
dilaton models. At high z the relative difference between 
both Hubble rates vanishes, as can be seen from Eq. CTl) 
as A —> 1 and £2 —> 0. 

We show in Fig.[2T]the deviation between the Einstein- 
frame and Jordan-frame density cosmological parame¬ 
ters. More precisely, for the Einstein frame we con¬ 
sider the cosmological “conserved” matter density p in¬ 
troduced in Ea. (IA13l) . Using Eo. (l37l) this is given by 


U m 


U m 

~T 


1 bn 

A(l-e 2 y~- 


(B2) 


paper, the background value of y < 10“ 3 is small, and 
the background value of K' ~ 1 implies that clusters are 
unscreened. We focus on processes that can happen in¬ 
side the hot gas of the clusters and consider the two-body 
scattering processes involving either electrons or nuclei. 
For temperatures of the gas less than T g < 10 keV, 
the particles are non-relativistic. We only consider K- 
mouflage functions K (y) whose small \ expansion starts 
at the cubic order. In this case, the three-point self¬ 
interaction of the scalar held is of order 


/• * ,,4 =(& p)(d r 5p) (d6p) 2 

Ci ~ M X - MR - 


(C2) 


where we consider a quasi-static background conhgu- 
ration. The energies of the two outgoing scalars are 
El ,2 ~ to/ whereas the spatial 3-momentum of the par¬ 
ticles is p r ~ JrrifTg. The scalar propagator gives a 
factor of 1 /to 2 in the non-relativistic limit and hnally we 
get that the scattering amplitude can be estimated as 


On the other hand, the radiation density parameter in 
the Einstein frame is f2( r ) = ff( r )/(l — e 2 ) 2 , from Eg. (1371) . 
and the Einstein-frame dark-energy density is then 




P 


m / 

MpiM 2 


X 3/2 V m f T a- 


(C3) 


^ ~ A — 1 ~ 

fJde = 1 f2(r) = ^ (B3) 

Again, we Hnd in Fig.(2l]that the differences between the 
Einstein-frame and Jordan-frame cosmological parame¬ 
ters are of the order of 1%. 


Appendix C: Unitarity constraints 

The K-mouflage models involve higher-order operators 
in the derivatives of <p and a coupling of the scalar held 
to matter /3. This coupling induces a Yukawa interaction 
of the scalar held with fermions, 

Cl = ( C1 ) 
Mpi 

where 5p = p — <p are the huctuations around a back¬ 
ground ip. Interaction terms in the scalar Lagrangian 
of the type M 4 x n in a background held configuration 
X imply the existence of the two-body scattering pro¬ 
cesses // —> pip at tree level, with the exchange of one 
scalar held ip. In quantum mechanics, unitarity of the 
scattering matrix requires that the scattering amplitude 
for this process should satisfy M < 167T. 

In the clusters that we consider in the main body of the 


For protons and neutrons at temperatures T g < 10 keV 
for j3 = 0.1 and using \ < 10 —3 , we hnd that < 

10 —2 , implying that unitarity is respected in the two- 
body processes. 

Terms of higher order in % can lead to processes such as 
// -A wp involving n scalars in the final state. The scat¬ 
tering cross section grows fast with the number of outgo¬ 
ing particles and can exceed the Froissart bound on the 
rr total "X: ppr- This violation is relaxed by classicalization 

0H] in the same fashion as in Galileon models ITst . 
where a classical lump sourced by the incoming energy of 
the two fermions in the center-of-mass frame is created. 
For K-mouhage models, this classical configuration has 
a typical size given by the K-mouhage radius Rk El 
and the scattering cross section becomes equal to the ge¬ 
ometrical cross section proportional to R 2 K . This process 
is analogous to the creation of black holes in high-energy 
collisions. 

K-mouhage models also satisfy a non-renormalization 
theorem analogous to the one for Galileons fH]. The 
quantum corrections going beyond the K-mouhage La¬ 
grangian are negligible when r > Aiy 1 , where Mk = 
K' x ! 4 M. Inside galaxy clusters, this is a short scale 
around 1 mm and we can completely neglect quantum 
corrections on cluster scales. 
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